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ABSTRACT 


This  note  presents  a  detailed  mathematical  analysis  of  a  multiple  - 
antenna  AMTI  radar  system  capable  of  detecting  moving  targets  over  a 
significantly  wider  velocity  range  than  is  achievable  with  a  single -antenna 
system.  The  general  system  configuration  and  signalling  strategy  is 
defined,  and  relationships  among  system  and  signalling  parameters  are 
investigated.  A  deterministic  model  for  the  target  return  and  a  statistical 
model  for  the  clutter  and  noise  returns  are  obtained,  and  an  optimum  pro¬ 
cessor  for  target  detection  is  derived.  A  performance  measure  applicable 
to  a  large  class  of  processors,  including  the  optimum  processor,  is  defined 
and  some  of  its  analytical  properties  investigated.  It  is  shown  that  an 
easily  implementable  sub-optimum  processor,  based  on  two-dimensional 
spectral  analysis,  performs  nearly  as  well  as  the  optimum  processor.  The 
resolution  and  ambiguity  properties  of  this  sub-optimum  processor  are 
studied  and  a  detailed  numerical  investigation  of  system  performance  is 
presented,  including  a  study  of  how  performance  varies  with  basic  system 
parameters  such  as  the  number  of  antennas. 
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A  THEORY  OF  MULTIPLE  ANTENNA  AMTI  RADAR 


1.  Introduction 

In  order  for  a  surveillance  radar  to  obtain  reasonable  resolution  in  both 
range  and  angle  and  also  have  good  coverage,  it  is  necessary  to  use  many  range  gates 
and  have  a  Doppler  processing  capability  for  each  range  bin.  Even  then,  the  wide 
spread  of  the  clutter  spectrum,  along  with  range  ambiguity  considerations  make  it 
extremely  difficult  for  a  high  speed  surveillance  airplane  to  detect  moving  targets 
embedded  in  ground  clutter.  This  latter  difficulty  may  be  greatly  alleviated  by  pro¬ 
cessing  signals  from  several  antennas  in  such  a  way  that  ground  clutter  returns  then 
are  essentially  cancelled.  Moving  targets  then  may  be  detected  even  if  their  Doppler 
frequencies  fall  within  the  main  beam  clutter  spectrum  of  the  return  received  at  each 
antenna. 

This  note  is  devoted  to  a  detailed  mathematical  exposition  of  the  theory  of  such 
a  multi-antenna  AMTI  radar.  The  major  goals  of  the  analysis  are  the  derivation  of 
an  optimum  processor  along  with  a  quantitative  evaluation  of  its  performance,  the 
formulation  of  an  easily  implementable  sub-optimum  processor  whose  performance 
is  close  to  that  of  the  optimum  processor  and  the  determination  of  a  good  choice  of 
the  system  parameters  such  as  number  of  antennas,  antenna  spacing,  etc. 

These  ends  will  be  accomplished  by  the  detailed  discussions  contained  in  the 
following  sections.  Section  2  contains  a  discussion  of  the  basic  system  parameters 
and  their  interrelations.  The  signalling  strategy  to  be  used  is  set  forth  and  an  ele¬ 
mentary  discussion  of  the  basic  clutter  cancellation  mechanism  to  be  exploited  by 
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the  system  is  included.  Section  3  continues  with  a  precise  description  of  the  target, 
clutter  and  receiver  noise  signals  that  are  received  by  the  radar.  The  desired  target 
return  is  described  in  a  deterministic  fashion  but  the  clutter  and  receiver  noise 
returns  are  modelled  as  random  effects  and  given  a  statistical  description  that  is 
suitable  to  the  forthcoming  analysis.  Section  4  addresses  the  problem  of  how  best 
to  extract  the  desired  signal  from  its  interfering  background  and  arrives  at  an  opti¬ 
mum  processor  for  this  purpose.  A  performance  measure  applicable  to  a  large 
class  of  processors  (including  the  optimum  processor)  is  defined  and  some  of  its 
analytical  properties  investigated.  Section  5  postulates  a  simple  sub-optimum  pro¬ 
cessor  and  gives  heuristic  as  well  as  mathematical  reasons  for  expecting  that  this 
processor  should  perform  almost  as  well  as  the  optimum  processor.  The  section 
concludes  with  a  short  discussion  of  the  resolution  and  ambiguity  properties  of  the 
sub-optimum  processor.  Section  6  concludes  the  body  of  the  note  with  a  detailed 
numerical  investigation  of  the  performance  of  the  optimum  and  sub-optimum  pro¬ 
cessors  and  provides  the  final  justification  for  adopting  the  latter  processor.  In 
addition,  an  investigation  of  various  system  configurations  leads  to  what  seems  to 
be  the  best  choice  of  the  basic  system  parameters.  The  major  results  obtained  are 
summarized  briefly  in  section  7  and  several  areas  in  which  more  work  should  be 
done  are  indicated. 


2 


2.  Geometry  and  Signalling  for  a  Multi-Antenna  AMTI  System 

a)  Basic  System  Configuration 

The  purpose  of  this  section  is  to  establish  constraints  among  some  of  the  basic 
geometric  and  signalling  parameters  of  an  AMTI  system  which  achieves  a  degree  of 
clutter  cancellation  by  processing  the  returns  from  several  antennas,  each  of  which 
has  collected  returns  independently  from  the  identical  clutter  area  on  the  ground.  The 
parameters  of  such  a  system  turn  out  to  be  so  closely  interrelated  that  (for  a  given 
aircraft  length  and  velocity)  a  choice  of  the  number  of  antennas  and  a  parameter 
relating  the  prf  to  the  aircraft  velocity  and  antenna  spacing,  essentially  determines 
the  remaining  system  parameters.  Hence  the  various  performance  criteria  (MTI 
gain,  unambiguous  range  interval,  velocity  resolution,  blind  speeds)  need  only  be 
investigated  for  variations  of  these  two  parameters. 

The  basic  geometry  for  the  problem  at  hand  is  depicted  in  Fig.  2. 1  where  the 
aircraft  is  assumed  to  be  moving  with  constant  velocity  v  in  the  positive  y  direction 
at  altitude  h.  A  linear  antenna  array  of  total  length  L  is  arranged  along  the  length 
of  the  aircraft.  By  appropriate  excitation  of  the  elements  of  the  array,  m,  essentially 
identical  synthetic  antenna  beams  are  formed  which  can  transmit  and  receive  independently 
at  different  times.  Each  such  antenna  beam  will  have  an  effective  aperture  length  denoted 
by  a  and  the  phase  center  will  be  equally  spaced  with  separation  d.  As  indicated  in 
Fig.  2.  2,  a  can,  if  desired,  be  made  larger  than  d  by  allowing  different  antennas  to 
share  some  of  the  elements  of  the  physical  array.  It  is  apparent  that 

L  =  (m-  1)  d  +  a  (2.  1) 

a  relationship  that  is  valid  whether  a  >  d  or  not.  The  forwardmost  phase  center  is 

th 

assumed  to  be  located  at  the  point  (0,0, h)  at  t=0which  means  that  the  i  phase  center 
(i=0,  ...  m-1)  is  located  at  point  f) ,  vt-id,  h)  at  time  t=0  where  d  denotes  the  separation 
between  successive  phase  centers. 


3 


11-6-13646 


L 


1  Fig.  2.2.  The  effective  antenna  apertures. 
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b)  Clutter  Cancellation  Principle 


It  is  quite  simple  to  obtain  a  quantitative  understanding  of  how  the  proposed 
system  cancels  ground  clutter  returns  without  at  the  same  time  cancelling  the  return 
from  a  moving  target.  The  basic  equation  for  this  purpose  relates  s.(t),  the  signal 
transmitted  on  the  i^1  antenna,  to  r.(t),  the  signal  received  on  the  i^  antenna  from 
a  stationary  point  scattering  target  on  the  ground.  Ignoring  amplitude  factors  which 
play  no  role  in  the  present  discussion  and  limiting  attention  to  the  case  m=2,  the 
equation  is 


Tj(t) 


S.  t  - 
1 


2Rt  (t) 
c 


1=0,  1 


(2.2) 


t  h 

where  c  denotes  the  velocity  of  light  and  R.(t)  denotes  the  distance  from  the  i 

antenna  to  the  ground  clutter  scatterer  (z  =0,  v  =  v  =  v  =0).  Reference  to  Fig.  2.  1 
°  o  ox  oy  oz 

shows  that  this  distance  is  given  by 


Rj(t)  = 


+  (y  +  id  -  Vt)2  +  h~ 
o  VJo 


11/2 


i= 0,  1 


(2.3) 


where  (xq,  yQ  0  )  is  the  location  of  the  stationary  ground  scatterer.  Now  if  the 
transmission  on  the  first  antenna  is  a  delayed  version  of  the  transmission  on 
the  zeroth  antenna,  s^(t)  =  s  (  t  -  A),  then 


r.(t)  =  s 


o 


2R.  (t) 


i=0  >  1 


(2.4) 


But,  it  is  easily  seen  that  if  A  =  d/v 
R.  (t  +  -)  =  R  ft) 

1  v  V  '  Ov  7 


(2.5) 


Kquation  (2.  2)  is  only  an  approximation  to  the  true  radar  return  but  it  can  be  shown 

that  this  approximation  is  valid  as  long  as  (— )^  <  <  X/R  ,  where  R  denotes  the 

c  max  max 

maximum  range  of  targets  of  interest.  This  condition  is  easily  met  for  most  cases 
of  interest. 
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and  it  follows  from  equation  (2.4)  that  the  signal  received  on  lagging  antenna  is 
identical  to  the  one  received  on  the  leading  antenna  except  for  a  known  delay. 

Since  this  holds  true  for  all  values  of  x^  and  y  ,  it  follows  that  the  returns  from 
all  fixed  clutter  scatterers  may  be  cancelled  by  delay  and  subtraction.  On  the 
other  hand  the  returns  from  a  mov ing  target  will  not  be  cancelled  by  this  delay 
and  subtraction  scheme  because  the  distance  the  transmitted  signal  must  travel 
from  the  leading  antenna  and  back  is  different  from  the  distance  the  transmitted 
signal  must  travel  from  the  lagging  antenna  and  back. 

The  technique  just  outlined  can  easily  be  extended  to  more  than  two  antennas 
and  can  be  applied  even  if  the  antennas  are  not  equally  spaced.  As  long  as  the 
time  offset  between  transmitted  signals  from  adjacent  antennas  is  equal  to  the 
time  necessary  to  fly  the  distance  between  the  antennas,  the  returns  at  each 
antenna,  due  to  fixed  clutter,  will  be  identical  except  for  a  delay.  This  mode 
of  operation  is  usually  referred  to  as  a  DPCA  mode  (Displaced  Phase  Center  Antenna)* 
and  this  nomenclature  will  be  used  hereafter  whenever  A,  the  time  between  successive 
transmissions  on  adjacent  antennas,  is  set  equal  to  d/v.  The  clutter  cancellation 
technique  just  described  is  basic  to  the  multi-antenna  AMTI  system  to  be  studied 
in  detail  in  the  remainder  of  this  report  and  a  clear  understanding  of  it  will  greatly 
facilitate  the  readers’  understanding  of  the  results  to  follow. 

c)  The  Illuminated  Range  Interval 

A  basic  requisite  of  the  clutter  cancellation  scheme  is  that  the  returns  from  the 

various  antennas  do  not  overlap  in  time,  so  that  they  can  be  sampled  independently 

without  crosstalk.  This  was  tacitly  assumed  in  the  discussion  above  which  used 
th 

the  fact  that  the  i  antenna  only  received  returns  that  were  produced  by  its  own 
transmission  and  that  these  returns  were  not  overlapped  by  returns  produced  by  trans¬ 
missions  from  the  other  antennas  at  other  times. 
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This  need  for  independence  of  antennas  imposes  some  additional  constraints 

on  the  signalling  strategy.  To  illustrate,  consider  the  signalling  scheme  in  Fig.  2.3, 

for  a  two  antenna  system,  where  s.(t)  represents  the  pulse  train  from  the  antenna, 

s^(t)  and  Sj(t)  are  assumed  identical  except  for  a  delay  of  d/v  and  1/r  is  the  prf  seen 

at  either  antenna.  Only  if  the  time  duration  of  the  return  from  a  single  pulse  is  less 

than  the  minimum  interval,  t  .  ,  between  the  transmitted  pulses  in  the  composite 

min  r  r 

pulse  train,  can  the  antenna  returns  be  sampled  independently.  This  limits  the 
maximum  range  interval  ^max  to  be  illuminated,  according  to  the  formula 
ct  . 


R 


max 


min 


(2.6) 


The  vertical  antenna  pattern  must  be  designed  so  that  only  this  range  interval  is 
illuminated. 

It  is  desirable  to  interleave  the  pulse  trains  so  that  R^ax  is  as  large  as  possible. 
Two  signalling  schemes  which  will  achieve  this  end  are  illustrated  in  Fig.  2.4.  In  each 
case,  the  pulse  trains  have  been  interleaved  so  that  in  the  steady  state  the  composite  pulse 
train  is  equispaced.  In  the  first  example, 


while  in  the  second  case 


Generalizing  from  these  two  examples,  it  is  seen  that  the  general  signalling  scheme 
satisfies, 


(2.7) 


(2.8) 


d 

v 


(2.9) 
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Fig.  2.  3.  A  two-antenna  signalling  scheme. 
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Fig.  2.  4.  Two  2-antenna  signalling  schemes. 
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where  N  is  the  number  of  pulses  transmitted  by  a  given  antenna  during  the  inter-antenna 

flight  time  d/v.  If  (2.  9)  is  imposed,  the  composite  pulse  train  will  be  equispaced  (in 

the  steady  state),  and  t  .  will  be  maximized  at  t/2  so  that 
J  min 

C  T 

R  =  £  75- 

max  2  2 


(2.  10) 


The  arguments  leading  to  (2.  9)  and  (2.  10)  easily  can  be  generalized  to  the  m-antenna 
case*  To  illustrate,  in  Figs.  2.5a  and  2.5b  two  possible  signalling  schemes  for  the 
4-antenna  case  have  been  indicated.  In  Fig.  2.5a,  N=1  where  N  has  the  same  meaning 
as  above,  and 


d 

v 


In  Fig.  2.5b,  N=2,  and 


7  =  <1+^>T 


In  both  cases,  the  illuminated  range  interval  is  given  by, 


(2.  11) 


(2.  12) 


R 


max 


c  r 
1  1 


The  generalization  of  (2.  12)  to  the  m-antenna  case  is  given  by 


7-[<n‘1>+^] 


(2.  13) 


(2.  14) 

The  considerations  that  led  to  the  constraintexpressed  by  equation  (2.  14)  can  be 
summarized  as  follows.  First,  it  was  required  that  the  transmitted  pulse  trains  from 
all  antennas  be  delayed  with  respect  to  each  other  in  such  a  way  that  all  antennas  transmit 
pulses  from  the  same  set  of  locations  in  space.  This  was  necessary  to  allow  for  the 
cancellation  of  fixed  clutter  returns.  Then, while  satisfying  this  clutter  cancellation 
constraint,  it  was  required  that  the  pulse  trains  be  interleaved  in  such  a  way  that  the 
range  swath  which  could  be  illuminated,  without  crosstalk  between  antennas,  was 
maximized.  The  resulting  illuminated  range  interval  is  given  by 
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Fig.  2.5.  Two  4-antenna  signalling  schemes. 
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(2.  15) 


R  -  . 

max  1  m 

By  designing  the  vertical  antenna  pattern  and  taking  advantage  of  the  signal  attenuation 
with  increasing  range  and  the  horizon  effect,  the  system  must  be  arranged  so  that 
appreciable  returns  are  received  only  from  this  portion  of  the  ground. 

One  further  point  should  be  noted.  It  is  true  that  crosstalk  between  antennas 
can  be  eliminated  by  limiting  the  illuminated  range  interval  as  indicated  above;  however, 
it  is  still  true  that  transmissions  from  a  given  antenna  can  be  received  by  the  other 
antennas.  The  only  thing  the  range  limitation  accomplished  was  to  insure  that  the 
returns  received  at  a  given  antenna  caused  by  transmissions  from  other  antennas 
do  not  overlap  in  time  with  the  returns  caused  by  the  given  antenna's  own  transmissions. 
These  non-selfgenerated  returns  were  ignored  in  the  above  discussion  of  the  clutter 
cancellation  principle  and  they  will  continue  to  be  ignored  in  the  remainder  of  the 
report.  It  is  clear  that  some  information  about  the  environment  the  radar  is 
probing  is  being  discarded  by  this  procedure  but  it  is  not  at  all  clear  that  this 
additional  information  can  be  used  to  increase  the  system's  clutter  cancellation 
effectiveness  appreciably.  This  is  a  point  that  should  be  the  subject  of  further 
investigation. 

d)  The  Sampling  Rate  Constraint 

It  is  convenient  at  this  time  to  introduce  a  spherical  polar -coordinate  system 
whose  origin  is  (0,0,  h).  This  coordinate  system  is  depicted  in  Fig.  2.  6  where  the 
origin  for  the  azimuth  angle  0Qis  chosen  to  be  the  (common)  antenna  bores ight  azimuth 
angle  G^  and  the  polar  angle  is  measured  from  the  negative  z-axis.  In  terms  of  these 
coordinates,  it  is  a  simple  exercise  to  see  that,  at  a  given  range,  R*  the  spread  of  Doppler 
frequencies  produced  by  ground  clutter  scatterers  (zq=0v  ^  =  v  =  vqz  =  0)  located  in 
the  main  lobe  of  any  of  the  m  antennas  is  approximately  given  by 
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Fig.  2.6,  Aircraft  centered  polar  coordinate  system. 
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W-X 


(2.  16) 


For  a  mechanically  scanned  antenna  whose  beamwidth  is  independent  of  0^,  W  would 
be  proportional  to  sin  0^.  But  for  the  electronically  scanned  antenna  under  consider¬ 
ation  the  beamwidth  varies  approximately  as  1/sin  0^,  and  hence  W  is  independent 
of  0^*  The  elevation  angle  o  to  be  used  in  equation  (2.  16)  is  the  angle  that  defines 
the  range  of  the  ground  clutter,  R’  cos  o'  -  h.  Now  it  is  well  known  that  any  radar 
using  repetitive  pulse  trains  will  cause  all  Doppler  frequencies  that  differ  by  the  prf 
to  appear  as  the  same  frequency.  If  the  prf  is  chosen  too  low  with  respect  to  the 
width  of  the  clutter  spectrum,  then  aliasing,  or  foldover  of  the  clutter  spectrum  will 
occur,  causing  the  effective  amount  of  clutter  that  must  be  combatted  at  a  given 
frequency  to  increase.  In  order  not  to  alias  the  clutter,  the  minimum  prf  of  each 
antenna  must  be  equal  to  the  spectral  width  of  the  clutter,  so  that  ~  -  W.  Actually 
to  avoid  aliasing  the  clutter  sidelobes,  the  minimum  prf  would  probably  be  chosen 
higher  than  W;  however,  the  prf  should  not  be  chosen  higher  than  necessary  because 

as  the  prf  is  raised,  the  maximum  illuminated  range  interval  R  decreases  as 
p  max 

per  equation  (2.  15).  The  dependence  of  the  prf  on  o’  might  appear  troublesome, 
since  the  implication  of  this  dependence,  together  with  (2.  1)  and  (2.  14)  is  that  the 
antenna  structure  would  have  to  change  as  a  function  of  this  angle.  However,  a  reason¬ 
able  procedure  would  be  to  design  the  system  by  imposing  the  constraint 


1 


2v 


-  =  a 
r  p  a 


(2.  17) 

where  g  >  1.  Equation  (2.  17)  guarantees  that  the  clutter  is  never  aliased  no  matter 


what  the  value  of  o’  is.  ^max  *s  now  not  always  the  largest  value  that  could  be  achieved 
but,  in  return,  spectral  sidelobes  of  the  clutter  have  been  decreased. 
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Another  disadvantage  of  aliasing  the  clutter  spectrum,  perhaps  more  impor¬ 
tant  than  the  effective  increase  of  clutter  power,  is  the  fact  that  stationary  ground 
targets  at  different  angles,  whose  Doppler  frequencies  differ  by  a  multiple  of  the 
prf,  become  indistinguishable.  Thus  the  capability  for  mapping  fixed  targets 
while  searching  for  moving  targets  is  lost.  If  the  prf  is  high  enough  to  prevent 
appreciable  aliasing  of  the  clutter  sidelobes  (i.e.  ^  2.5),  the  Doppler  filtering 
of  the  returns  received  at  any  one  of  the  antennas  would  produce  a  synthetic  aperture 
ground  map  whose  cross-range  resolution  is  determined  by  the  physical  beam  width 
of  an  individual  antenna  divided  by  the  number  of  returns  processed.  It  will  turn 
out  (see  section  5)  that  the  signal  processing  to  be  employed  for  AMTI  includes 
Doppler  filtering  of  the  returns  received  at  each  antenna  by  means  of  a  discrete 
Fourier  transform.  Thus  for  a  high  enough  prf  the  AMTI  processing  includes,  as 
an  intermediate  step,  essentially  all  the  processing  necessary  for  forming  a 
synthetic  aperture  ground  map. 

If  the  system’s  clutter  rejection  capability  were  so  great  than  the  additional 
clutter  power  caused  by  aliasing  were  of  no  consequence,  and  if  ground  mapping 
were  of  no  concern,  consideration  could  be  given  to  dropping  the  prf  requirement 
(2.  17).  The  prf  could  then  be  lowered,  with  a  consequent  increase  in  the  illuminated 
range  interval.  However,  as  indicated  in  an  example  below,  it  appears  that  quite 
sufficient  range  coverage  can  be  obtained  even  while  the  prf  is  high  enough  to  pre¬ 
vent  aliasing.  Finally,  even  if  the  prf  requirements  (2.  17)  is  dropped,  a  lower 
bound  on  the  prf,  and  therefore  an  upper  bound  on  the  illuminated  range  interval, 
is  imposed  by  the  DPCA  constraint  (2.  14). 
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In  the  ensuing  discussion  and  examples  it  will  be  assumed  that  clutter  is  not 
aliased,  and  the  prf  is  chosen  at  least  equal  to  (and  perhaps  higher  than)  the  clutter 
spectral  width. 

Effects  of  Constraints 

The  constraints  that  have  been  established,  now  can  be  brought  together  and 
solved  for  their  joint  effect.  The  parameters  that  have  been  of  concern  so  far 
include: 

d  =  the  inter-antenna  spacing 

a  =  the  length  of  each  individual  antenna  aperture 
L  =  the  overall  length  of  the  antenna  array 
v  =  aircraft  speed 
1/t  =  the  prf  at  each  antenna 
Rmax  “  the  unambiguous  range  interval  of  the  system 
m  =  the  number  of  antennas 

N  =  the  number  of  pulses  transmitted  by  a  given  antenna 
before  the  next  antenna  in  line  moves  to  the  original 
spatial  position  of  that  given  antenna. 
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The  constraints  between  these  parameters  have  been  stated  in  (2.  1),  (2.  14) 
(2.  15),  and  (2.  17).  These  constraints  may  be  summarized  as: 


L  =  (m-1)  d  +a  (geometry)  (2.  18) 

vr  =  d/(N  -  ~~)  (clutter  cancellation)  (2.  19) 

vr  =  a/ (2/i)  (sampling  rate  to  prevent  aliasing)  (2.  20) 

R^ax  =  2  m  (unam^8uous  range  interval)  (2.21) 


If  we  take  L,  v,  and  [i  as  given,  then  for  a  specified  m  and  N,  these  equations 
can  be  solved  explicitly  for  the  remaining  parameters.  The  results  are: 


,  (N-l)  m  +  1 _ 

C  (N-l)  (m-T)  m  +  (1  +  2/j)  m-1 


(2.22) 


2m M 

(N-I)  (m-l)ln  +  (1  +^[m-'l 


(2.23) 


T 


m _ 

(N-l)  (m  - 1 )  m  +  JT ~ 2[i)  m  - 1 


L 

v 


(2.24) 


1  CL 

^max  (N-l)  (m-1)  m  +~(1~+  2 fu)  m-1  ~2v 

It  can  be  observed  from  (2.25)  that  ^max  *s  largest  for  m=2,  N=l,  and  decreases 
monotonically  with  increasing  m  or  N.  This  statement  imples  that  unless  other  considerations 
(such  as  MT1  gain  or  velocity  resolution)  dictate  otherwise,  a  two-antenna  system,  with  a 
simple  alternate  pulsing  of  the  antenna  would  be  chosen.  It  will  be  shown  later  that  after 
consideration  of  these  factors,  m=2,  N=1  remains  a  reasonble  choice  of  system  para¬ 
meters. 

Examples  of  Parameter  Values 

This  section  will  be  closed  with  a  few  examples  of  the  parameter  values  dictated 
by  (2.22)  -  (2.25).  It  is  assumed  in  the  examples  that  ^-l. 
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First  leave  N  arbitrary,  and  consider  the  two  antenna  case.  The  results 


are: 


(2.  2d) 


a  TN+3  L 


(2.27) 


(2. 28) 


R 


1  cL 
max  Tn+3  Tv 


(2.29) 


Note  that  R 


max 


is  a  monotonically  decreasing  function  of  N,  and  is  maximum  when  N=F 


=  I  cJk 

max  5  Tv 


(2.30) 


The  inter -antenna  spacing  is  1/4  of  the  aperture  size,  implying  a  need  for  antennas 


which  overlap  and  must  share  elements.  For  a  Mach  1  (1100  ft/sec. )  aircraft  with  a 


10  foot  array,  the  illuminated  range  interval  is  about  170  miles. 

Form  (2.  26)  and  (2.  27)  it  is  seen  that  the  ratio  d/a  of  spacing  to  individual 

antenna  size,  is  less  than  1  (requiring  overlapped  antennas  unless  N  :>  3.  However,  for 

large  N,  the  R  decreases. 

&  max 

For  a  comparison  to  the  two-antenna  case,  consider  the  example  of  m=4,  N=l. 
The  results  are 


=  cL 
max  11  Tv 


(2.31) 


The  inter-antenna  spacing  is  1/8  of  the  length  of  each  antenna.  For  a  given  L/v, 


the  illuminated  range  interval  is  slightly  less  than  half  of  that  in  the  corresponding  two-antenna 
case. 
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3.  Target  and  Clutter  Returns 


The  purpose  of  this  section  is  to  construct  a  mathematical  model  of  the  signals 

received  by  an  m-antenna  airborne  array  used  in  the  manner  described  in  the  previous 

section.  Three  types  of  received  signals  will  be  considered:  1)  The  signal  produced 

by  a  desired  target  which  is  defined  to  be  a  moving  point  object  located  on  or  above 

the  ground,  2)  The  signal  produced  by  non-moving  ground  clutter  and  3)  The  signal 

produced  by  receiver  front-end  noise.  The  basic  structure  of  the  radar  receiver  to 

be  used  is  sketched  in  Fig.  3.  1.  The  returns  produced  by  each  transmitted  pulse 
th  th 

from  the  i  antenna  are  received  by  the  i  antenna,  corrupted  by  front-end  noise 

and  then  passed  to  a  filter  that  is  matched  to  the  (common)  subpulse  shape.  The 

output  of  the  matched  filter  is  then  range  gated  with  each  gate  producing  a  sample 

of  complex  video  from  the  same  range  for  each  pulse  transmitted.  The  samples 

obtained  in  this  manner,  along  with  the  samples  from  the  other  antennas,  are  then 

passed  to  a  processor  whose  nature  will  be  specified  later  and  a  decision  as  to  the 

presence  or  absence  of  a  desired  target  in  each  range  gate  is  made. 

In  the  sequel,  attention  will  be  focused  on  a  single  range  gate.  The  samples 

obtained  from  this  gate  when  n  pulses  are  transmitted  on  m  antennas  can  be  arranged 

as  un  nm  dimensional  column  vector  denoted  by  r.  More  precisely,  let  the  complex 

th  th 

sample  received  on  the  i  antenna,  i=o,  . . .  m-1,  due  to  the  k  pulse,  k=o, . . .  n-1, 
transmitted  by  that  antenna  be  denoted  by  r.^.  Then  r  is  a  column  vector  whose 
components  are  arranged  in  the  order  r^,  r^  ,  ...  r^  r^,  rn9  •••  n-l’### 


rm-l,0’  rm-l,  1*  '  "  rm-l,n-l* 

The  received  vector  r  can  be  represented  in  the  form, 


(3.1) 
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COMPLEX  RECEIVED 
SAMPLES 


1 1-1-15151 


DECISION 


NOISE 


RECEIVED  SAMPLES 
FROM  OTHER  ANTENNAS 


Fig.  3.1.  The  basic  radar  receiver. 
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where  c  and  n  denote  the  received  vectors  produced  by  ground  clutter  and  receiver 
noise  respectively.  The  vector  s  denotes  the  signal  received  from  the  desired 
moving  target  normalized  in  such  a  way  as  to  bring  out  explicitly  the  dependence 
of  this  return  signal  on  the  angular  location  of  the  target  in  the  antenna  beam,  the 
radar  cross  section  of  the  target,  and  the  transmitted  energy.  Accordingly,  G(6q) 
denotes  the  two-way  voltage  gain  of  each  of  the  antennas,  evaluated  at  the  azimuth 
0Q  (measured  from  the  antenna  boresight  azimuth  0^  as  per  Fig.  2.  6)  of  the  desired 
target.  It  has  been  assumed,  as  is  usually  the  case,  that  the  vertical  beam  pattern 
is  sufficiently  broad  so  that  the  antenna  gain  is  essentially  independent  of  the  target’s 
elevation  angle  relative  to  the  boresight  elevation  angle  for  all  targets  in  the  illuminated 
range  interval.  The  scale  factor  Aq  is  a  constant  that  accounts  for  the  radar  cross 
section  of  the  desired  target  and  the  signal  attenuation  due  to  range.  The  factor  E 
represents  the  transmitted  energy  per  pulse.  In  the  following  subsections,  detailed 
descriptions  of  the  received  vectors  £,  c  and  n  will  be  given. 

a)  The  Desired  Signal  Vector  £ 

The  purpose  of  this  sub-section  is  to  compute  the  components  of  the  vector  s. 
These  components  are  the  received  samples  produced  when  the  transmitted  signals 
are  returned  by  an  arbitrarily  located,  moving  point  scatterer. 

Assume  that  the  location  of  the  scattering  object  at  t=o  is  the  point  (xo,y^,z^) 

and  that  the  rectangular  coordinates  of  its  velocity  are  (vqx,  v  ,  VQZ)*  It  follows, 

as  shown  in  Fig.  2.  1,  that  the  location  of  the  scattering  object  at  time  t  is  (x^  +vox  t, 

th 

y^  +  vQy  t ,  zq  +  vQzt).  The  complex  signal  r.(t)  received  on  the  i  antenna  is  given 
by  equation  3.  2, 
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(3.2) 


where  s.(t)  denotes  the  complex  envelope  of  the  signal  transmitted  by  the  i  antenna 


normalized  to  unit  average  energy  per  pulse  transmitted  and  R.(t)  denotes  the  distance 

th  r 

between  the  i  antenna  and  the  moving  object.  The  amplitude  factor JEA^G  (0q)  has 


been  deleted  in  equation  (3.  2)  because  it  already  has  been  taken  into  account  by  equation 
(3.  1).  Reference  to  Fig.  2.  1  shows  that 


(3.3) 


This  last  expression  is  much  too  complicated  to  be  of  use  in  the  following  analysis; 


however,  useful  approximations  to  equation  (3.3)  can  be  obtained  by  expanding  it  in  a 
double  power  series  about  t=0,  d=o  and  retaining  only  the  first  few  terms.  These 
approximations  turn  out  to  be  valid  as  long  as  the  observation  time  is  sufficiently  small 
and  the  range  to  the  target  is  large  compared  to  d.  When  only  the  constant  and  linear 
terms  of  the  power  series  are  retained,  the  resulting  approximation  to  R.(t)  is  given  by 


o 


(Tx  v 

Mo  ox 


4-  y  (v  -  v)  +  (z  -  h)  v 
*0  ov  O  oz 


1 1  +  y  id 

;  J  o  . 


(3.4) 


o 


o 


where 


(3.5) 


and  vq  denotes  the  target’s  radial  velocity  at  t=o. 
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Making  use  of  the  polar  coordinates  defined  in  Fig.  2.6,  equation  (3.4)  can 


be  rewritten  in  the  form, 

R . (t)  a?  R  +[v  "  v  sin  cos  +  O  t  +  id  sin  o  cos  (0  +  0  ) 

1  0  [  0  0  boj  0  bo 

th 

The  signal  transmitted  by  the  i  antenna  is  assumed  to  be  a  uniformly  spaced 

coherent  pulse  train  consisting  of  n  pulses.  In  symbols,  the  signal  transmitted  on  the 
.th  . 

1  antenna  is 


where  u(t)  denotes  the  complex  envelope  of  the  transmitted  subpulse  (usually  either  a 

simple  pulse  or  a  linear  fm  pulse-compression  waveform)  normalized  so  that 

*  2  th 

u  (t)  dt  =  1  and  a.^  denotes  the  complex  amplitude  of  the  k  pulse  transmitted 

from  the  antenna.  The  transmitted  signal  s.(t)  is  normalized  to  unit  average  energy 

per  pulse  by  requiring 


1  = 


nm 


m-  1 
2 

i=o 


s.(t) 


dt  = 


nm 


m-  1 
2 

i=o 


n-  1 
2 

k=o 


ik 


As  before,  t  denotes  the  prf  seen  at  one  antenna  and  A  denotes  the  transmission  delay 
between  antennas.  The  transmission  delay  A  will  always  be  set  equal  to  d/v  to  achieve 
DPCA  operation  as  discussed  in  the  previous  section.  For  the  time  being,  however,  it 
will  simplify  notation  to  allow  A  to  be  arbitrary.  It  is  further  assumed  that  the  duration 
of  u(t)  is  smaller  than  either  A  or  r  so  that  no  pulses  in  the  train  overlap  in  time. 
Combining  equations  (3.2),  (3.  6),  and  (3.7)  yields  the  following  expression  for  the  signal 
received  by  the  1*1  antenna, 


(3.6) 


(3.7) 


(3.8) 
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2R 

r,  «)  =  exp  [j2»  £  (t  -  — exp  [-J2-T  |  (vo  -  vao)  t] 

n  ~  1  2R  r  ?  v  1  r  *] 

Z  a.k  u  (t  -  iA  -  kr  -  — ^-)  exp  - j2rr  (id  -^)  exp  |  -j2 n  £  (iA  +  kr) 
k=o  L  J  L  J 

where  the  target’s  apparent  velocity  is  defined  by, 

v  =  v  sin  o  cos  (9,  +  0  ) 
ao  o  v  b  o' 

In  arriving  at  equation  (3.  9)  it  has  been  assumed  that  only  the  gross  antenna -to -target 

2R 

delay  o/c  has  an  appreciable  effect  on  the  subpulse  envelope  u(t). 

The  output  of  each  of  the  airborne  antennas  is  next  passed  through  a  filter  matched 
to  the  subpulse  shape  and  then  range  sampled  before  being  processed  further.  It  is 
possible  to  justify  this  initial  subpulse  matched  filtering  and  range  sampling  operation 
on  theoretical  grounds;  however,  the  argument  is  fairly  involved  mathematically  and  does 
not  give  any  real  insight  into  the  operation  of  a  multi-antenna  airborne  MTI  system.  For 
this  reason,  as  well  as  the  fact  that  subpulse  matched  filtering  seems  a  natural  thing 
to  do,  this  point  will  not  be  pursued  further  and  the  formal  discussion  of  the  processor 
will  begin  with  the  range  samples  derived  in  the  above  manner. 

ttl 

The  output  of  the  subpulse  matched  filter  connected  to  the  i  antenna  is  given  by 

yi(t)  =  J*  r.  (n)  u*  (77  -  t)  exp  [-j2rr  S.  -  t)]  dr? 

2RT 

This  output  is  next  sampled  at  the  times  — — +  iA  +k r,  yielding  the  desired 
complex  samples 


(3.9) 


(3.  10) 


(3.  11) 
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n- 1 


(3. 12) 


2r* 

s •!  =  y.  (—7—  +  iA  +  kr),  i=o,  . . .  m-1,  k=0, . . . 
ik  1  c 

This  yields  nm  samples  corresponding  to  the  target  returns  from  a  single  range  R'. 

Samples  from  other  range  gates  are  obtained  by  incrementing  R’  by  an  integer  multiple 
of  the  range  resolution  of  a  subpulse.  After  some  tedious  but  straightforward  manipulation, 
equation  (3.  12)  can  be  combined  with  equations  (3.  9)  and  (3.  11)  to  yield  the  expression, 


Sik 


at*  expWo>  [-i2"  X  [kT  <v„  -  *ao>  +  ivao  £  -  i)  +  *  Vu  j 


(3. 13) 


y  (R  -  R\  v  -  v  ) 

'll  o  o  ao 

where  (^,77),  the  ambiguity  function  of  the  subpulse,  is  defined  by 

yu  =  v"u  ^  u*  (4  0  exP  (-^  X  W)  dt  (3-  14) 

and  where  denotes  a  carrier  phase  angle  that  is  independent  of  i  and  k.  In  arriving 

at  equation  (3.  13),  it  has  been  assumed  that  the  scatterer’s  range  is  less  than  Rmax* 

Equation  (3.  13)  can  be  simplified  in  several  ways.  First,  the  requirement 
for  DPCA  operation  (A  =  d/v)  causes  the  second  term  in  the  exponential  to  vanish. 

Next,  the  fact  that  in  any  practical  situation,  the  Doppler  frequencies  involved  are 
much  smaller  than  the  reciprocal  of  the  subpulse  duration  makes  the  approximation 
(R()  -  R\  vQ  -  v^o)  (Rq  -  R’,  o)  a  reasonable  one.  These  simplifications 
result  in  the  final  equation, 


sik~  aik  exP0%)exP  [-i2"  l  [*r  (vQ  -  vao)]+  i  d^2. 

Finally,  it  will  be  assumed  that  the  desired  target  is  in  the  range  resolution  cell 
being  sampled,  |R  -  R'  |<  l/Wu  where  is  the  bandwidth  of  a  subpulse.  It  then 
follows  that  (Rq  -  R\  o)  a:  ^(o,  o)  =  1. 


(3.  15) 
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b)  The  Clutter  Vector  c 

The  ground  clutter  will  be  modelled  as  a  continuum  of  stationary  scatterers. 
th 

The  signal  received  on  the  i  antenna  from  a  small  patch  of  such  scatterers  c.(t),  could 
be  expressed  mathematically  by  setting  vq  =  o  in  equation  (3.  9)  and  multiplying  it  by 
a  suitable  amplitude  factor;  the  result  is, 

c  .(t)  =^E  G(8)  A  (R,  0)  Rd  Rd  0  exp  j 2tt  £  (t  -  -^)J  exp  £j2rr  — t~j 

^  3ik  u(t-iA)-kr-fl)  exp[-j2rr  |  (id  -^-)|exp[  +  kx)]  (3.  16) 

where  now, 

v  =  v  sin  o  cos  (6,  +  0) 
a  o 

and  R  and  8  denote  the  range  and  azimuth  of  the  scattering  patch.  The  elevation  angle 
of  the  patch  is  related  to  R  by  means  of  the  relationship  R  cos  co  ~  h.  The  function 
A  (R,  0)  is  the  clutter  amplitude  desnity  per  unit  area.  As  mentioned  before,  this 
amplitude  factor  takes  into  account  the  attenuation  with  range  as  well  as  the  clutter 
cross  section 

Unfortunately,  equation  (3.  16)  neglects  certain  effects  which  are  vital  to 
an  understanding  of  the  clutter  cancellation  capabilities  of  the  multiple -antenna  AMTI 
system  under  consideration.  As  mentioned  in  the  previous  section,  the  key  phenomenon 
the  processor  to  be  described  exploits  is  the  fact  that,  in  the  DPCA  mode,  each  antenna 
receives  essentially  the  same  clutter  signal  from  the  ground.  The  key  word  here  is 
’’essentially".  If  all  the  antenna  gains  were  truly  identical  and  if  the  pulses  transmitted 
on  different  antennas  were  exactly  the  same,  it  turns  out  that  the  processor  would 
ahieve  perfect  clutter  cancellation  and  system  performance  would  be  limited  by  receiver 
front-end  noise  alone.  Of  course,  this  is  never  the  case  in  practice  and  so  it  becomes 

important  to  model  in  some  way  the  effects  of  small  differences  in  the  clutter  returns 
received  by  the  different  antennas. 
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This  end  will  be  accomplished  by  assuming  that  the  gain  of  the  i^  antenna 
is  given  by  G.(0),  and  that  the  transmitted  subpulses  on  the  i^  antenna  are  given 
by  u.(t).  These  functions  are  assumed  to  differ  from  their  nominal  values  G($  and 
u(t)  only  in  "small"  ways.  This  notion  of  smallness  will  be  made  more  precise 
later.  It  now  follows,  that  the  clutter  return  from  a  small  patch  of  scatterers 
located  at  (R,  0)  is  given  by  equation  (3.  16)  modified  by  replacing  G(®  by  G.(0)  and 
u(t)  by  u.(t).  The  total  clutter  return  is  obtained  by  integrating  the  so  modified 
equation  (3.  16)  over  all  R  and  0.  The  resulting  signal  is  then  passed  through  a 


2r* 

filter  matched  to  the  nominal  subpulse  shape  u(t)  and  sampled  at  the  times  — —  4*  iA  +  kr 


to  yield  the  clutter  samples  c.^_.  The  result  of  this  somewhat  tedious  calculation  is. 


c . 


ik 


E  aik  J’  jA(R,0)  exp 


tMR.0) 


G.(0)  Vi  (R  -  R',o) 


2v  T 

exp  (j2rr  — ^ -  k)  RdRdB 


where  (R,  E)  is  the  carrier  phase  angle  corresponding  to  in  equation  (3.  13)  and 
v-  (s  7)  is  the  cross -ambiguity  function  defined  by, 


^  (5,^)  =  JV(t)  u*  (t  C)  exp 


dt 


-j2r  x  rjt 

All  the  simplifications  leading  to  equation  (3.  15)  except  the  one  requiring  the  scatterer 
to  be  in  the  range  resolution  cell  being  sampled  have  been  made  in  arriving  at  equation 


(3.  18) 


(3.  19) 


(3. 18). 


The  clutter  returns  c  given  by  equation  (3.  18)  are  critically  dependent 

ik 

on  the  detailed  way  in  which  A(R,  0)  exp  j>^(R,  £)  varies  with  R  and  0;  in  fact,  this 
function  usually  is  quite  noise-like  in  character.  For  this  reason,  it  is  common  practice 
to  model  the  clutter  as  noise  and  to  use  a  statistical  description  of  the  clutter  return 
to  design  and  evaluate  systems  for  combatting  such  clutter.  One  such  statistical  para¬ 
meter  is  the  clutter  covariance  matrix  which  will  play  a  vital  role  in  the  subsequent 
analysis  and  which  can  be  calculated  from  equation  (3.  18)  once  some  assumptions  as  to 
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the  statistical  behavior  of  A(R,  0)  exp  j\J/(R,  0)  are  made.  Two  common  and  very 

2  3 

useful  assumptions  are  as  follows:  ’ 

1.  A(R,  0)  and  \T/(R,  0)  are  statistically  independent  and  >F(R,  0)  is  uniformly 
distributed  in  the  interval  (-tt,tt).  From  this  it  follows  that  the  mean  of  A(R,  0)  exp  j\J(R,  0) 
is  zero,* 


£  |a(R ,  e)  exp  j  >£(R ,  0)J  =  0 

2.  Scatterers  located  at  different  locations  in  space  are  statistically  independent. 
Because  of  equation  (3.  20),  this  assumption  implies  the  equation, 

A(R,  0)  A  (R\  0')  exp  j  |V(R,  0)  -  \E'(R' >  0')j 

=  cr  (R,  0)  6  (R-R’)  6  (0-0’) 


(3.  20) 


(3.21) 


where  <7  (R»  0)  denotes  the  average  clutter  ’’cross-section”  density  per  unit  area  and 
S(x)  denotes  the  Dirac  delta  function. 

Given  these  assumptions,  it  is  now  an  easy  matter  to  calculate  the  clutter 

kk’ 

covariance  matrix  C  whose  elements  c..f  are  defined  by 

n  J 

kk'  r  *  1 

cii’ =  eLcikc  ivj 

and  the  result  is  given  by 

2  V  7“ 

c^’=E  aika*.,k,  f  Ja(R,0)  F.  (R,0)  F.,(R,  0)  exp[j2n  ^  -  (k-k’)]  RdRd0 

where 

Fi  (R,  0)  =C.(®y.  (R-R’,o) 


(3. 22) 


(3. 23) 

(3.  24) 


*The  symbol  £  denotes  the  expectation  operator. 

■**■*00113  functions  are  used  here  and  later  in  the  description  of  receiver  noise 
strictly  for  mathematical  simplicity.  This  artifice  gives  correct  results  as 
long  as  the  actual  correlation  lengths  or  times  are  small  compared  with  the 
resolution  of  the  system. 
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The  clutter  matrix  C  is  to  be  formed  out  of  its  elements  c.^  according  to  the  scheme 
given  by, 


C  = 


00 

01 

0,  m- 1 

c 

C 

C 

11 

12  •' 

1 ,  m-1 

(3.25) 

Hermitian 


m-l,n-l 


where  the  submatrices  C..t  are  defined  by, 

n  J 


C  f  = 
ii 


00 

c.., 

n 


"  iV 


11 
C  iif 


0,  n-  1 

c 

n 


1,  n- 1 

c 


Hermitian 


n- 1  ,n- 1 

c 

ii 


(3.  26) 


c)  The  Noise  Vector  n 

The  last  effect  that  must  be  accounted  for  is  that  due  to  receiver  front-end 

th 

noise.  Denoting  the  noise  voltage  introduced  by  the  i  antenna  by  n.(t),  it  is  seen  from 
equations  (3.  11)  and  (3.  12)  that  the  received  noise  samples  n  ^  are  given  by, 


n.  =  n.  IP)  u  (r?  -  iA  -  kT  - 
lk  J  i  v  7 


2R’ 


)  drj 


(3.27) 


A  statistical  description  of  this  noise  can  be  given  by  making  the  usual  assumption  that 
the  real  receiver  noise  is  a  zero-mean,  white  Gaussian  process  with  a  power  spectral 
density  of  Nq  watts/Hz.  This  corresponds  to  a  complex  Gaussian  noise  with  a  auto" 
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correlation  given  by, 


n.(t)  n*  (f) 


=  4  N  6  (t-t’) 
o  ' 


(d.28) 

In  addition  it  will  be  assumed  that  the  noises  introduced  by  different  front-ends  are 
statistically  independent  which  means  that 


nj(t)  n?,  (t') 


=  4  N  6. .,6  (t-t’) 

O  11  v  ' 


(3.29) 


where  <5..,  denotes  the  Kronecker  delta.  It  now  follows  easily  that  the  covariance  matrix 


of  the  noise  samples  has  elements  given  by, 
P 


ik  11  i’k' 


=  4  N  <5..,  <5,  ,  t 
o  u*  kk’ 


(3.30) 


which  means  that  tne  receiver  noise  covariance  matrix  is  of  the  form  4  N  I  where  I  is 

o 

the  (nm)  x  (nm)  identity  matrix.  Finally,  assuming  that  the  receiver  noise  is  statistically 
independent  of  the  clutter,  it  follows  that  the  combined  covariance  matrix  tor  clutter 
plus  noise,  A  is  given  by, 

A=C  +  4N  I  (3.31) 

n  o 

The  matrix  will  be  referred  to  as  the  interference  matrix  in  the  sequel, 
d)  Summary 

It  has  been  shown  that  the  received  signal  vector  given  by  equation  (3.  1) 
consists  of  a  desired  signal  vector  s  whose  components  are  given  by  equation  (3.  15), 
a  clutter  vector  whose  components  are  given  by  equation  (3. 18)  and  a  receiver  noise 
vector  whose  components  are  given  by  equation  (3.27).  The  clutter  vector  has  been 
described  statistically  by  means  of  its  covariance  matrix  C  whose  elements  are  given 
by  equation  (3.23),  and  the  receiver  noise  covariance  matrix  is  equal  to  4NqI.  The 
covariance  matrix  of  the  combined  (clutter  plus  noise)  interference  is  given  by  eq.  (3.31) 
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4. 


Data  Processing 


Now  that  a  detailed  description  of  the  data  collected  by  the  radar  has  been 
given,  it  is  possible  to  address  the  question  of  now  best  to  process  these  data  to 
separate  a  moving  target  from  its  background  of  ground  clutter  and  receiver  noise. 

The  present  section  will  argue  that,  for  several  reasons,  a  logical  choice  for  this 
processor  is  a  linear  filter  followed  by  an  envelope  detector.  A  performance 
measure  for  such  processors  will  be  defined  and  the  processor  that  maximizes 
this  performance  measure  will  be  derived.  This  optimum  processor,  as  well  as 
its  associated  performance  measure,  turns  out  to  depend  on  the  covariance  matrix 
of  the  interference.  Since  this  matrix,  as  derived  in  the  previous  section,  is  quite 
difficult  to  use  for  numerical  evaluations,  certain  additional  assumptions  will  be 
made  in  order  to  simplify  its  structure.  In  addition  to  simplifying  numerical  work, 
these  assumptions  will  prove  to  give  a  great  deal  of  insight  into  the  nature  of  the 
interference  and  the  structure  of  the  optimum  processor.  The  next  task  of  this 
section  will  be  to  define  the  "improvement  factor”  for  an  arbitrary  linear  processor. 
This  factor  will  be  normalized  in  such  a  way  as  to  make  it  quite  simple  to  compare 
AMTI  systems  having  different  parameters  such  as  number  of  antennas  or  number 
of  pulses.  It  will  be  shown  that,  subject  to  some  additional  assumptions,  the  improve¬ 
ment  factor  factors  into  two  terms,  one  associated  with  temporal  (pulse-to-pulse  on 
a  given  antenna)  processing  and  the  other  associated  with  spatial  (antenna-to-antenna) 
processing.  This  result  is  valuable  because  it  better  defines  the  structure  of  the 
processor  and  because  it  gives  a  quantitative  measure  of  the  relative  contribution 
spatial  and  temporal  processing  make  to  the  total  improvement  factor.  Finally,  it 
will  be  demonstrated  that  the  optimum  processor  can  usually  be  considered  as  a 
cascade  of  a  purely  temporal  processor  and  a  purely  spatial  processor. 
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a)  The  Optimum  Linear  Processor 

One  of  the  most  well  justified  and  commonly  used  techniques  for  defining 

the  structure  of  the  data  processor  to  be  used  for  extracting  a  desired  signal  from 

4  5 

a  background  of  interference  is  statistical  decision  theory.  ’  This  theory  requires 
that  a  complete  statistical  description  of  the  interference,  in  this  case  the  ground 
clutter  and  receiver  noise,  be  available.  Because  both  clutter  and  receiver  noise 
can  be  thought  of  as  the  superposition  of  a  large  number  of  small,  statistically 
independent  voltages,  it  is  common  practice  to  invoke  the  central  limit  theorem 
and  assume  that  the  receiver  noise  and  clutter  voltage  are  Gaussian  and,  therefore, 
completely  specified  in  a  statistical  sense  by  a  knowledge  of  their  means  and  covariances. 
Given  this  assumption,  it  is  well  known  that  statistical  decision  theory  as  applied  to 
the  problem  of  extracting  a  signal  that  is  known  except  for  its  complex  amplitude  from 
background  interference  leads  to  a  data  processor  that  consists  of  a  suitably  chosen 
linear  weighting  of  the  received  data  followed  by  an  envelope  detector.  The  output  of 
the  detector  is  compared  with  a  threshold  and  ’’target  present”  is  announced  when 
the  thresnold  is  exceeded. 

This  result  can  be  stated  in  precise  mathematical  terms  by  means  of  the 

notation  r,  s  and  A  defined  in  the  previous  section.  Using  this  notation,  the  processor 

- 1 

discussed  above  can  be  shown  to  calculate  the  quantity  f  =  |  r  A  £  |  and  compare 

it  with  a  threshold.  Since  the  signal  vector  s  depends  on  the  unknown  target  velocity 

vq  and  target  azimuth  and  elevation  (0^,^)  through  the  quantity  v^  =  v  sin  oq  cos  (^+0^* 

the  detection  statistic  i  must  be  calculated  for  all  values  of  v  and  v  that  are  of  interest 

o  ao 

and  a  separate  threshold  comparison  must  be  made  for  each  pair  of  these  values.  In 
other  words,  each  threshold  comparison  tests  the  hypothesis  and  there  is  a  target 
present  whose  parameters  vq  and  v  are  equal  to  those  values  used  in  computing  J2. 

*The  asterisk  denotes  the  conjugate  transpose  of  a  matrix. 
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The  processor  just  described  is  a  special  case  of  a  general  linear  processor 
which  is  defined  to  be  a  device  that  computes  £  r*  w  where  w  is  an  arbitrary 
nm -dimensional  vector,  and  compares  it  with  a  threshold.  The  weights  w  =  A  1  s_ 
which  can  be  derived  via  statistical  decision  theory  under  the  assumption  that  both 
receiver  noise  and  clutter  are  Gaussian  actually  have  an  important  optimality 
property  which  does  not  require  the  assumption  that  the  underlying  interference  is 
Gaussian.  In  order  to  state  this  property,  the  signal -to -interference  ratio  (SIR)  of 
the  linear  processor  is  defined  to  be 

2 

(output  for  target  alone  present  ) 

output  variance  for  no  signal  present 


EA2  I  G  <e0)  I2 


#  2 
Is  wf 

w*  A  w 


(4.1) 


Now  it  is  a  simple  exercise  in  applying  the  Schwarz  inequality  to  show  that  the  SIR 
defined  by  equation  (4.  1)  is  maximized  when  w  =  A  s  and  that  this  maximum  value 
is  given  by 

p  =  EA  2  I  G  (e  )  I2  S*  A'1  s 
hmax  o  1  '  o  1  -  - 


(4.2) 


Thus,  the  linear  processor  whose  weighting  vector  w  is  given  by  w  =  A  1  s  is  optimum 
in  the  sense  that  it  maximizes  the  SIR  given  by  equation  (4.  1) 

In  the  light  of  the  above  discussion  it  seems  reasonable  to  use  a  linear  processor 
with  the  optimum  weights  given  by  w  =  A  ^  s  to  separate  the  desired  signal  from  its 
background  interference.  This  is  the  course  that  will  be  followed  in  the  sequel  even 
though  it  shortly  will  be  necessary  to  introduce  a  number  of  additional  statistical 
effects  which  presumably  may  be  non-Gaussian  into  the  description  of  the  clutter. 
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The  linear  receiver  will  now  no  longer  necessarily  be  the  one  dictated  by  statistical 
decision  theory  but,  it  will  still  maximize  the  SIR  over  the  set  of  all  linear  receivers 
as  long  as  the  weights  are  chosen  according  to  the  formula  w  =  A  ^  where  A  is  the 
covariance  matrix  of  the  interference  obtained  by  averaging  over  all  statistical 
effects,  Gaussian  and  non-Gaussian. 

Before  discussing  these  additional  statistical  effects, it  is  important  to  note 
that  because  a  linear  processor  is  to  be  used,  there  is  no  reason  to  transmit  anything 
but  equal  amplitude  zero-phase  pulse  trains.  In  other  words,  no  loss  of  generality 
is  incurred  by  assuming  a.^  =  1  in  the  remaining  discussion.  This  conclusion  comes 
about  because  of  the  fact  that  the  receiver  weights  w.^  and  the  transmitted  weights  a.^ 
always  appear  as  a  product  w  a.^  in  the  equation  for  the  SIR  (4.  1).  Thus  it  follows 
that  anything  that  can  be  accomplished  by  weighting  on  transmit  can  equally  well  be 
accomplished  by  weighting  on  receive.  Since  the  latter  is  easier  to  accomplish  than 
weighting  on  transmit  it  will  henceforth  be  assumed  that  a.^  =  1. 

The  motivation  for  the  introduction  of  additional  statistical  effects  lies  in  the 
complexity  of  the  clutter  covariance  matrix  C  which  is  needed  to  compute  the  optimum 
weights.  Reference  to  equations  (3.  22)  and  (3.  23)  reveals  that  the  elements  of  C 
depend  on  the  cross-ambiguity  functions  x.  (£»  0)  and  the  gain  functions  G.(0),  i=0, .  .  . 
m-1.  So  far  in  the  discussion  it  has  been  tacitly  assumed  that  these  quantities  were 
known.  This  is  not  really  the  case  in  practice  because,  for  example,  the  y.  (£,0)  only 
differ  from  one  another  because  of  small,  presumably  uncontrollable,  departures  of 
the  transmitted  subpulses  u.  (t)  from  the  nominal  pulse  u  (t).  Similarly,  even  if 
all  the  antennas  were  designed  to  have  the  same  nominal  gain  G(9),  slight  fabrication 
errors  can  result  in  departures  of  the  individual  gains  from  the  nominal  gain.  These 
effects  often  can  be  neglected  but  the  present  system's  clutter  rejection  properties 
depend  heavily  on  having  these  quantities  equal  and  thus  the  effect  of  small  departures 
from  equality  must  be  taken  into  account. 
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A  convenient  way  to  do  tnis  and  one  that  leads  to  quite  tractable  final  results 
is  to  assume  that  the  departures  of  the  y.  (R  -  R\0)  and  the  G.(0)  from  their  nominal 
values  are  statistical  in  nature  and  then  to  make  certain,  hopefully  reasonable, 
assumptions  about  these  statistics.  The  simplest,  non-trivial  assumptions  that  can 
be  made  concern  the  first  few  moments  of  the  pertinent  statistics.  To  this  end,  the 
following  assumption  will  be  made  about  the  second  moment  of  F.  (R,  0), 


F.(R,0)F*  (R ,  0) 


=  (1  I  f (r , e) 


(4.3) 


where  S'..,  =1  -  6..,  and  F(R,  0)  =  G  (0)  \u  (R  “  R\0). 

In  addition,  it  will  be  assumed  that  F.(R,  0)  is  statistically  independent  of  the  random 
scattering  mechanism  producing  the  clutter.  The  quantity  O'  appearing  in  equation  (4.3) 
is  presumed  to  have  a  magnitude  that  is  small  compared  with  unity.  It  is  possible, 
of  course,  to  make  more  involved  and  perhaps  more  realistic  statistical  assumptions 
about  F.(R,  0)  than  those  just  made  but  to  do  so  would  to  some  extent  bring  back  the 
original  dependence  of  C  on  the  details  of  F.(R,0)  which  the  statistical  assumptions 
were  supposed  to  remove.  The  efficacy  ot  assumption  (4.  3)  in  predicting  actual 
system  performance  will  require  validation  by  experiment  or  simulation. 

The  clutter  covariance  matrix,  equation  (3.23)  must  be  recomputed  now  in  order 
to  take  account  the  effect  ot  the  additional  statistical  parameters  that  have  been  introduced. 
Because  of  the  assumed  statistical  independence  of  F.(R,0)  and  the  ground  scattering 
mechanism,  the  existing  equation  (3.  23)  can  be  interpreted  as  given  the  conditional 
covariances  of  the  clutter  given  F^,  0);  all  that  remains  to  be  done  is  to  average 
over  the  remaining  random  variables.  With  reference  to  equation  (4.3),  the  result 
of  this  operation  easily  is  seen  to  be, 


cf,  =E  <1  -a6..IxfJ,a  (R,0)  |  F  (R,  0) 


expjj2n 


2v  r 

-x-  <k  -  k’> 


RdRde 


(4.4) 
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This  definition  of  a  implied  by  equation  (4.3)  is  somewhat  unsatisfactory 
because  it  does  not  give  any  insight  into  how  to  estimate  a  in  any  given  situation. 
This  can  be  remedied  to  some  extent  by  noting  that  for  i  ^  if, 

£E  J  Ja  (R,0)  I  F.(R,0)  -  F  (R,e)  | 2  RdRdG 


=  2E  J  Ja  (R,0)  I  F  (R,0)  |2  RdRd0 
-  2  (l-»)  E  J  Ja  (R,9)  |  F  (R,  9)  | 2  RdRdG 

-  2E  cy  J  J  cr  (R.  9)  |  F  (R,  G)  |2  RdRdG  (4.5) 

so  that,  2 

J  c  J  Ja  (R,0)  I  Ft  (R,0)  -  F..  (R,0)  |  RdRdG 

J  Jct  (R,6)  |  F(R,6)  | 2  RdRdG  (4.6) 

Equation  (4.  6)  can  be  used  to  derive  a  more  useful  estimate  for  a  by  first  assuming 
that  ct  (R,  0)  is  a  constant  and  then  noting  that, 

eJ'J  |  Xf  (R-R')[g.(0)  -Gi,(G)]+G.,(G)[xi  (R-R')  -  ^  (R-R')]  | 2  RdRdG 
J  J  |  Xu  (R-R')  |2  I  G  (6)  I2  RdRdG 


eJJ  I  xu<R-R,> 


G.(0)  -  G.,(e)]  +  G(G) 


Xj  (R-R’)  -  Xj.(R-R') 


J*  J  |  Xu  (R-R’)  I2  |  G(6)  I2  RdRdG 


RdRdG 


ej  I  G.  (G)  -G.,(G)  | 2  d g  e  J  I  Xj  (R-R’)  -yr  (R-R’)  |2RdR 

J  I  G(Q)  |2  dG  J  |  xu  (R-R1)  I2  RdR 

where  second  order  terms  such  as  j^G.(0)  -  G(G)JJg.(0)  -  G.,(G)  J  have  been  ignored, 
G.  and  y.  have  been  assumed  statistically  independent  and  it  has  been  assumed  that 

c[Gj  (6)]  =  G(6)  and  (R-R’)]=XU  (R-R')- 
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Equation  (4.7)  provides  a  means  of  estimating  a  in  terms  of  average  percentage  antenna 
errors  and  average  percentage  pulse  mismatch  errors. 


Equation  (4.4)  can  be  written  in  a  simpler  form  by  the  introduction  of  some 
new  notation.  To  this  end,  the  input  clutter  power  per  pulse  is  defined  by 

pc  =  c*f  =  E  J  Ja  (R,8)  I  F  (R,0)  |2RdRd6 
and  the  temporal  correlation  function  H(x)  is  defined  by 


2  r  2v  T  1 
H(x)  =  E  P J 1  j'J’cKR.e)  |  F(R,  6)  I  exp  [2jTT  -±-  XJ  RdRd0 


In  terms  of  this  notation,  (4.4)  can  be  written  as 


_  IrVr  '  — 

c  ...  =  P  (1  -  a  6.., )  H  (k  -  k') 

11  c  '  ii’ '  '  ’ 

Finally,  the  entire  interference  covariance  matrix  A  -  C  +  4N  I  will  be  written  in 

o 

terms  of  the  normalized  matrix  M  defined  by 


A»PCM 

It  is  easily  seen  that  the  elements  of  M  are  given  by 


kk’  _ 
m  =  p 
u  r 


nc 


6.., 

n 


6kk,  +(1  -  fffi..,)  H  (k  -  k’) 


where 


P 


nc 


4N 

o 

P 


c 

is  defined  to  be  the  input  noise -to -clutter  ratio  per  pulse. 


(4.8) 


(4.9) 


(4. 10) 


(4.11) 


(4. 12) 
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c)  The  Improvement  Factor 

The  output  signal -to- interference  ratio  p  defined  by  equation  (4.  1)  obviously 
can  be  considered  as  a  figure  of  merit  for  a  given  AMTI  system  and  systems  having 
different  parameters  (number  ofpulses,  number  of  antennas)  and/or  processors  could 
be  compared  by  comparing  their*  SIR's.  Such  comparisons  can  be  simplified  somewhat 
by  the  introduction  of  the  notion  of  an  improvement  factor  which  measures  the  SIR 
gain  that  is  obtained  from  a  given  system. 

The  improvement  factor  can  be  defined  by  first  noting  that  the  SIR  for  an 
arbitrary  linear  processor  given  by  equation  (4.  1)  can  be  written  in  the  form, 

EAo  i  c,(eo)  i2  i  £*  y  r 

p  p  * 

c  w  M  w 

where  use  has  been  made  of  equation  (4.  11).  Next,  it  is  observed  that  the  quantity 

9  i  2  -  1 

EA^  I  F(Q^)  |  can  be  identified  as  an  input  signal -to -cl utter  ratio  per  pulse, 


(4.13) 


P: 


in 


G(0O) 


(4.  14) 


The  improvement  factor  I(vq,  v  )is  now  defined  to  be  the  ratio  p/p.^  and  is  given 


bv, 


I(v  ,v  ) 
'  o  ao 


*  ,  2 
s  w 


w  M  w 


(4.  15) 


where  the  dependence  of  the  improvement  factor  on  the  desired  target’s  velocity  and 
location  parameters  v^  and  v^q  (through  the  components  of  s)  has  been  indicated  explicitly 
in  the  notation.  Several  points  should  be  noted  about  this  definition  of  improvement 
factor.  First,  although  p  was  defined  as  an  input  signal  to  clutter  ratio  and  not 
an  input  signal -to- interference  ratio,  these  two  quantities  are  almost  equal  in 
any  practical  situation  where  usually  the  input  clutter  power  far  exceeds  the  input 
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noise  power.  The  next  point  to  bear  in  mind  is  that  two  systems  should  always  be 


compared  with  the  desired  target  at  the  same  azimuth  for  otherwise,  p.^  will  be 
different  for  the  two  systems  and  the  improvement  factor  will  not  give  a  fair  comparison. 
In  the  same  vein,  equation  (4.  5)  shows  that  systems  having  different  antenna  apertures, 
a,  will  have  different  p.^'s  because  their  antenna  gains  and  beamwidths  will  be  different. 
This  last  effect  turns  out  not  to  be  very  important,  however,  because  the  system 
constraints  discussed  in  section  2  (see  equation  (2.  23))  do  not  allow  a  to  vary  significantly 
for  N=1  which  turns  out  to  be  the  most  interesting  case.  Finally,  it  should  be  noted 
that  the  improvement  factor  for  the  optimum  linear  processor,  I mSLX  (vn»  vaQ)  can  ^ 
written  in  the  form, 


I  (v  ,  v  )  =  s*  M  *  s  (4.  lb) 

max'  o  ao7  -  - 

An  extensive  numerical  investigation  of  the  improvement  factor  will  be  presented 
in  a  later  section. 

d)  Spatial -Temporal  Factorization 

The  remainder  of  this  section  will  be  devoted  to  showing  that  as  long  as  the  clutter 
power  is  much  larger  than  the  receiver  noise  power,  the  improvement  factor  for  a 
large  class  of  interesting  linear  processors  factors  into  a  product  of  two  terms,  one  of 
which  gives  the  improvement  due  to  temporal  (pulse-to -pulse)  processing  and  the  other  of  which 
gives  the  improvement  due  to  spatial  (antenna-to -antenna)  processing.  This  is  an 
important  result  because  it  facilitates  the  assessment  of  the  relative  efficacy  of 
increasing  the  number  of  pulses  or  the  number  of  antennas.  In  addition,  it  will  be 
shown  that  subject  to  the  same  assumption  the  optimum  processor  has  the  property 
that  it  factors  into  two  independent  processors  one  performing  purely  temporal 
processing  and  the  other  performing  purely  spatial  processing. 
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The  establishment  of  these  factorizations  is  greatly  aided  by  the  introduction 
of  the  Kroneker  product  of  two  matrices.  The  Kroneker  product^A  ®  B  of  an  n  x  m 
dimensional  matrix  A  with  elements  a.^  with  a  pxq  dimensional  matrix  B  with  elements 
b^s  is  defined  by  the  equation 


a 


11 


a12B 


H..B 

ik 


A  ®  B  = 


a21B 


a2kB 


a  ,B  ...  a  B 

ml  mn 


(4.  17) 


Note  that  the  matrix  A  &  B  is  np  x  nq  dimensional.  The  sequel  will  make  use  of 
several  easily  verified  identities  involving  Kroneker  products.  They  are, 


(A  ®  R)~  1  =  A"1  ®  B~ 
if  A  and  B  are  non- singular , 


(4.  18) 


(A  8  B)  (C  8  D)  =  (AC)  ®  (BD) 
provided  that  the  products  AC  and  BD  are  defined, 

(A  S  B)*  =  A*  ®  B* 


(4.19) 


(4.20) 


Inspection  of  equation  (3.  15)  now  shows  that  the  vector  £  can  be  written  in 

+ 

the  form  s  =  s  ^  where  and  respectively  denoting  the  temporal  and  spatial  parts 
of  s,  are  given  by, 


+  The  unimportant  carrier  phase  angle  has  been  dropped  from  this  equation. 
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(4.21) 


s* 

-s 


r  2 (v  -  v  )t  -i  r  2(v  -  v  )  -il 

1,  exp  j^-j2rr  - —  ^  — - J  ,  exp  Jj-j2rr - 2_ — (n-l)|j 

1,  exp  (-j2n  -2_),  ...  exp  jjj^rr  Q  (m-l)jj 


In  addition,  inspection  of  equation  (4.i0)shows  that  the  matrix  P  *  C  can  be  written 
in  the  form  P  *  C  =  M,  ®  M  where  the  elements  of  M  are  given  by 

C  L  S  t 


Mkk'  =  H  (k  "  k') 


k,k'  =  0,  1 


.n-1 


and  the  elements  of  Mg  are  given  by 

m|.(  =  (1  -  Of  6  ,  i,  i*  =  0, 1,  . .  .  m-1 


It  will  now  be  assumed  that  the  clutter  power  is  large  enough  so  that  receiver 
noise  effects  can  be  ignored.  The  precise  mathematical  definition  of  what  is  meant 
by  this  assumption  is  that  be  small  compared  to  the  eigenvalues  of  P  *C.  This 
assumption  implies  that  C.  The  factorization  of  the  imorovement  factor  now 

can  be  established  for  all  separable  processors;i.e. ,  processors  whose  weights 
satisfy, 

w  =  wA  ®  w 

-  -t  -  s 

where  w^  is  an  arbitrary  1  x  n  vector  and  Wg  is  an  arbitrary  1  x  m  vector.  Given 
these  assumptions  and  making  use  of  the  identities  (4.  20;  and  (4.21),  the  improvement 
factor  can  be  expressed  as, 

2 

|  (s  f  ®  s  )*  (w  ®  w  )  | 

I(v  ,  V  )  =  - - - - - - - =5 - 

°  ao  (w '  ®w  )*(M  ®M  )(w,®w  ) 

t  S  t  S  “  t  —  s 

I  (s*  ®  s*)(wt®  Wg)  I  2 
(w*  ®  w£ )  f(Mtwt)®  (Mgw  ^ 


(4.2  2) 


(4.  23) 


(4.24) 


(4.  25) 
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2 


I  £*  ”t  1 2  I  -*  -si 

<“*  Mt  ^  <”*  Ms  V 


=  I,  (v  ,  V  )  I  (v  ) 

t  o  ao  s  o 

where  the  temporal  and  spatial  improvement  factors  are  given  by 


I  St*  wt  I2 

’t^o-  vao>  =  TT - 

St  Mt^t 


and 


I  (v  )  =  — - 
s  o  * 


*  ,  2 

s  w 

s  — s  1 


s  Ms 

-  s  s  —  s 


(4.26) 


(4.27) 


(4.28) 


Equation  (4.  26)  exhibits  the  desired  spatial -temporal  factorization  of  the  improvement 
factor. 


It  is  easy  to  see,  using  the  Schwarz  inequality  that  the  maximum  achievable 

improvement  for  a  separable  linear  processor  is  given  by, 

I  (v  ,  v  )  =  (s*  M,  1  s, )  (s*  M  1  s  ) 
max'  o  ao7  -t  t  -t7  -s  s  -s' 

and  the  associated  weighting  vectors  are, 


st 


w  =  M  1  s 
— s  s  —  s 


(4.29) 


(4.30) 


The  notation  I  (v  ,  v  )  used  in  equation  (4.  29)  tactitly  assumes  that  the 
max  v  o  ao  M  J 

maximum  achievable  improvement  with  a  separable  processor,  subject  to  the  assumption 
M  ^  P  *  C  is  comparable  to  the  improvement  achievable  with  the  optimum  processor,  w 
=  M  1  s,  subject  to  the  same  approximation.  This  follows  from  the  fact  that  the 
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optimum  processor  is  approximately  separable  as  may  be  seen  from  the  following 
equations, 

m'1  s-(g)"1  s  =  (Mt®  Ms)_1  (st®  ss) 
c 

=  (Mt‘1®M‘1)(st«Ss) 

=  K'  !*> 

where  use  has  been  made  of  identity  (4.  18) 


(4.31) 
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5. 


The  Double  Discrete  Fourier  Transform  (DDFT)  Processor 


This  section  will  be  devoted  mainly  to  the  derivation  of  an  approximation 

to  the  optimum  processor  discussed  in  the  last  section.  As  will  be  shown,  the 

7 

approximate  processor  is  a  device  that  first  performs  m  n-point  DFT’s  on  the 
pulses  received  from  each  antenna.  The  resulting  data  is  then  processed  in  an 
antenna-to -antenna  fashion  by  means  of  n  m-point  DFT's.  The  magnitude  of 
each  of  the  nm  complex  numbers  derived  in  this  manner  is  then  compared  with 
a  threshold  and  ’’target  present”  is  announced  for  each  threshold  crossing.  The 
DDFT  frequency  resolution  cell  for  which  the  threshold  crossing  occurred  will 
be  shown  to  give  a  measure  of  the  target's  radial  velocity  vq  and  its  apparent 
velocity  v^Q  =  v  sin  oq  cos  (0  +  0q)  which  gives  some  information  about  and 
0  .  The  usefulness  of  the  approximate  processor,  which  will  be  shown  to  perform 
within  a  few  db  of  the  optimum  processor,  lies  in  the  fact  that  a  DFT  can  be  performed 
with  extreme  ease  and  rapidity  on  a  digital  computer  by  means  of  the  well  known 
fast  Fourier  transform  (FFT)  algorithm.  The  optimum  processor,  on  the  other 
hand,  requires  much  more  complicated  and  time  consuming  computation.  The 
section  will  continue  with  the  derivation  of  an  approximate  formula  for  the  improve¬ 
ment  factor  of  the  DDFT  processor  and  conclude  with  a  discussion  of  its  spatial  and 
temporal  resolution  and  ambiguity  properties. 


a)  Derivation 


The  first  step  in  deriving  the  desired  approximate  processor  is  to  introduce 


the  derived  data  d.^,  i=0,  . . .  m-1;  k=0, . .  .n-1  which  is  obtained  from  the  original 
data  by  means  of  a  DDFT.+ 


l  m-1  n-1 


2  2 
p-0  q=0 


(5.1) 


s  t 

+  Equation  (5. 11)  is  actually  a  "windowed"  DDFT  because  of  the  weights  w  and  w  but  this  fact 
will  not  be  reflected  in  the  nomenclature. 
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s  t 

where  w^  and  are  respectively  arbitrary,  non-zero  spatial  and  temporal  weighting 

factors.  Since  the  DDFT  is  well  known  to  be  an  invertible  operation  (the  r.^  can  be 

recovered  from  the  d.^.  by  a  formula  that  is  identical  to  equation  (5.1)  except  that  the 

sign  of  the  exponents  is  reversed),  it  follows  that  the  optimum  linear  processor  for 

extracting  the  desired  signal  from  the  derived  data  d.^  must  perform  just  as  well  as 

the  optimum  linear  processor  based  on  the  original  data  r  The  processor  based 

on  the  derived  data  is  defined  by  the  expression  SL  =  |  d*  D  *  v  |  where  d  is  a  vector 

composed  of  the  derived  samples  d.^  in  the  same  manner  that  r  was  composed  of  the 

samples  r.^.  D  denotes  the  covariance  matrix  of  the  derived  data  whose  elements 

d^t  are  defined  by 
11  J 


=  fi  (d d*  ) 

x  lk  1  k  ' 


(5.2) 


when  clutter  and  receiver  noise  alone  are  present  in  d.^.  The  vector  v  is  the  DDFT  of 
the  desired  signal  vector  and  has  elements  v.^  defined  by 


m- 1 


ik 


p=0 


n-1 

2 

q=0 


pq 


wS  w*  exp  (j 2tt  2L  )  exp  (j2rr  3^) 
pq  m  n 


(5.3) 


It  should  be  noted  that  the  DDFT  defined  by  equation  (5.  1)  can  be  implemented 
as  a  cascade  of  a  temporal  and  a  spatial  DFT  in  any  order  desired.  Furthermore,  for 
the  case  of  only  two  antennas,  the  spatial  DFT  merely  consists  of  adding  and  subtracting 
n  pairs  of  complex  numbers. 

The  next  task  will  be  to  demonstrate  that  the  matrix  D  is  approximately  a  diagonal 
matrix.  An  intuitive  argument  for  why  this  should  be  so  can  be  constructed  along 
the  following  lines.  The  temporal  part  of  the  DDFT  takes  a  set  of  n  pulses  from  each 
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antenna  and  frequency  analyzes  them  into  n  frequency  resolution  cells.  If  this 
frequency  analysis  were  ideal,  the  only  contribution  to  a  given  cell  would  come 

from  those  clutter  scatterers  whose  apparent  velocities  produced  Doppler 
frequencies  that  lay  within  the  bounds  of  that  frequency  resolution  cell.  It  now 

follows  from  the  initial  assumption  (equation  3.21)  that  clutter  scatterers  located 
at  different  points  on  the  ground  scatter  independently,  that  the  outputs  of  different 
frequency  resolution  cells  -  because  they  are  looking  at  different  frequency  intervals, 
and,  therefore,  different  regions  on  the  ground  -  must  be  independent  also.  A  DFT 
does  not  perform  an  ideal  frequency  analysis  as  required  by  the  above  argument  because 
the  filters  it  uses  have  sidelobes.  Nevertheless,  a  properly  "windowed”  DFT  can  have 
filtering  properties  that  come  reasonably  close  to  ideal  which  implies  that  the  above 
conclusion  is  true  to  within  a  good  approximation.  The  argument  so  far  has  established 
that  all  elements  of  D  that  refer  to  the  correlation  between  different  temporal  frequency 
cells  are  approximately  zero. 

The  next  step  in  the  argument  is  to  note  that  the  output  of  any  of  the  n  spatial 
DFT's  is  significant  only  in  that  spatial  resolution  cell  corresponding  to  zero  frequency. 
This  follows  from  the  fact  that  the  use  of  the  DPCA  mode  makes  the  clutter  highly 
correlated  from  antenna -to -antenna  and,  therefore,  essentially  a  dc  voltage  that 
shows  up  primarily  in  the  zero-frequency  resolution  cell.  This  last  argument 
establishes  the  fact  that  any  term  in  D  which  refers  to  the  correlation  between  different 
spatial  resolution  cells  must  also  be  approximately  zero.  Combining  these  two  conclusions 
leads  to  the  final  result  that  D  is  approximately  diagonal. 

The  elements  of  D  can  be  evaluated  explicitly  by  reference  to  equation  (4.  6), 

(4.  12),  (5.  1)  and  (5.2)  with  the  result 
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p 

c 

nm 


m-1 

2 

p,p'=0 


n-1 

2 

q,qf=0 


qq'  s  s*  t  t* 

m w  w  ,  w  w  f 
pp.  P  P  q  q 


exp  O^Ei^expO^3^) 


(5.4) 


which,  after  some  straightforward  algebra,  can  be  rewritten  in  the  form 


r>  S  t  -  p  .S  -t 

=  P  p  e..,e  +Ef..,  f,  ,, 

c  rnc  li  Kk  li  kk 


(5.5) 


where 


s  1 

e..,  = — 
li  m 


m-1 
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p=0 


| 2  exp  p  1  1 
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(5.6) 
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k-k’ 
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|2  exp  (j2rrp  ilL) 

# 
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(1-cy) 
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2 

_p=0 

w®  exp  (j2n  Ei) 

'm-1 

2 

Lp=o 

w^  exp  (j2n  Ei) 
p  m 

(5.8) 

fkk-  I  ^’6)  I2fi  (l 


2v  T 
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2v  T 
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(5.9) 
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and  where, 


fiw = i 


n_1  t 

2  w  exp  (j2rrqx) 
q=0  q 


1 


n- 1 


f9  =  ^  2 

2  n  q=0 


w*  |  2  exp  (j2nqx) 


(5. 10) 


It  now  will  be  shown  that,  as  long  as  the  spatial  and  temporal  weights  are  properly 
chosen,  the  covariance  matrix  D  is  approximately  a  diagonal  matrix.  In  this  connection 
it  is  important  to  note  that  the  performance  of  the  optimum  processor  based  on  the  derived 
data  is  independent  of  the  choice  of  temporal  and  spatial  weights  as  long  as  none  of  these 

g 

weights  are  zero.  To  this  end,  a  uniform  set  of  spatial  weights  w^  will  be  chosen  with 
the  result  that, 


e..t  =6.., 

li  if 


(5.11) 


and 

[“+»>  6io  siojaii- 

The  temporal  weights  are  to  be  chosen  in  such  a  manner  that  the  functions 

f(jX)  and  f2(x)  defined  by  equation  (5. 10)  behave  roughly  as  indicated  in  Fig.  (5.  1) 

That  this  sort  of  behavior  can  be  achieved  by  a  suitable  choice  of  weights  w^  is  well 

Q 

known  from  the  theory  of  window  functions  used  in  spectral  analysis.  For  example, 
if  uniform  weights  are  used  both  f^(x)  and  f£(x)  are  the  f°rm  sin  nx/sin  x  which  exhibits 
the  desired  behavior  but  has  rather  large  sidelobes.  The  theory  of  window  functions 
shows  how  to  depart  from  uniform  weights  so  as  to  achieve  low  sidelobes  without 
appreciable  broadening  of  the  main  lobe. 


(5.  12) 


48 


18-6*13652 


Fig.  5.1.  Desired  behavior  of  fj(x)  and  f2(x). 
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Assuming  that  f^(x)  and  f9(x)  have  negligible  sidelobes  and  assuming  that  the 

1  ^  1  £  2 

weights  have  been  normalized  so  that  -2  I  w  =  1  it  now  follows  that 

n  q=0  q 


In  addition,  it  is  obvious  that  the  integral  in  equation  (5.  9)  is  essentially  zero 
unless  k=k'  which  means  that  fj^,  is  approximately  given  by 

ft  _  ft  . 
kk'  kk  kk' 

fkk  =  n  J  MR’e>  I  F  (R-e)  I2  I  fi  <|  +  ~T~)  I2  RdRde 

The  next  fact  to  be  established  is  that  the  desired  signal  vector  v  has  only 
one  non -negligible  component.  Intuitively,  this  result  is  reasonably  obvious  because 
equation  (3. 15)  shows  that  the  desired  signal  is  a  two-dimensional  complex  exponential 
which  will  tend  to  be  almost  orthogonal  to  all  the  exponential  factors  in  the  DDFT 
except  the  one  whose  frequencies  match  those  of  the  desired  signal  most  closely.  This 
result  can  be  established  more  precisely  by  substitution  of  equation  (3.  15)  into  equation 
(5.  3)  with  the  result, 

^  K  WQ  eXP[J2nql|-^  (Vo'Vao| 

1  ™-‘  Pw*  M  VI 

r  U  expr*- 

=  fnm  f  |  ["—  -  (v  -v  )1-  2  exp  [j2TTp  (-L  -  ^  )] 

y  1  |_n  \  '  o  ao  J  m  „  r  |_J  vm  A.  v  'J 


(5.  13) 


(5.  14) 

(5.  15) 


(5.  16) 
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Because  of  the  assumed  behavior  of  f^(x),  the  first  term  of  equation  (5.  16)  is  approximately 


equal  to  zero  for  all  values  of  k,  k=0 


l  •  •  • 


closest  to  an  integer.  Similarly,  since  the  second  term  of  equation  (5. 16)  is  of  the  form 

sin  mx/sin  x,  it  is  approximately  equal  to  zero  for  all  values  of  i,  i=0,  . . .  m-1  except 

2d  v 

that  value  which  makes—  -  — -  closest  to  an  integer.  It  follows  that  the  vector  v, 

m  A  v  - 

as  stated,  has  only  one  non -negligible  component. 

It  now  follows  from  the  above  discussion  that  the  optimum  weighting  vector 
D  *  v  has  only  one  component  whose  value  differs  appreciably  from  zero.  This  means 
that  the  optimum  processor  defined  by  SL  =  |  d*  D  *  v  |  has  an  extremely  simple  structure 
in  that,  for  each  hypothesis  as  to  the  value  of  the  target’s  unknown  velocity  vq  and  apparent 
velocity  v^q  to  be  tested,  only  that  value  d.^.  corresponding  to  the  ’’closest”  spatial -temporal 
resolution  cell  need  be  compared  to  a  threshold.  The  ’’closest"  resolution  cell  corresponding 

to  a  given  vq  and  v^  is  defined  by  those  indices  i’  and  k’  that  make  the  quantities 

i  2d  Vo  ,  k  2t  x  ,  .  , 

—  -  __  —  and  _  -  (v  )  closest  to  integers, 

m  A  v  n  A  o  ao  b 

This  argument  establishes  the  assertion  that  the  DDFT  processor  is  an  approx¬ 
imation  to  the  optimum  processor.  The  goodness  of  this  approximation  will  be  established 
by  means  of  numerical  examples  in  a  later  section  where  it  will  be  seen  that  the  approximate 
processor  performs  essentially  as  well  as  the  optimum  processor  in  all  cases  of  interest. 


b)  The  Improvement  Factor 


The  task  of  this  section  is  to  derive  a  simple,  approximate  formula  for  the 

improvement  factor  I  (v  ,  v  )  of  the  DDFT  processor.  The  function  I  (v  ,  v  ) 
r  max'  o  ao  F  '  o  ao 

was  defined  to  be  equal  to  the  ratio  of  output  SIR  to  input  signal -to-clutter  ratio  and 


max'  o 


is  easily  seen  to  be  given  by, 


(5.17) 


c 
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Making  use  of  the  same  approximations  used  in  deriving  the  DDFT  processor,  it 
is  seen  that  equation  (5. 17)  is  approximately  equal  to 


.  v  )  - 


Kv 


<v  a0"  d.^.T'/p 

11  c 


(5. 18) 


where  i'  and  k'  denote  the  indices  of  the  spatial  and  temporal  resolution  cells  "closest" 

to  v  and  var).  Equation  (5. 18)  will  be  evaluated  explicitly  by  means  of  suitable 
l  i2  k'k'  k'k' 

approximations  for  |  v.,.,  |  and  d .  The  variance  d  is  given  by  equation  (5.5) 

1  iC  11  11 

and  reference  to  equations  (5. 11)  through  (5. 15)  shows  that  the  only  quantity  that  requires 
approximation  is  the  integral  given  by  equation  (5. 15).  This  integral  will  be  evaluated 
approximately  by  means  of  the  following  considerations. 

The  width  of  a  range  resolution  cell  is  assumed  to  be  small  enough  so  that  &  in  the 
integral  is  approximately  equal  to  the  constant  p* ,  where  d  is  the  elevation  angle  to  the 

center  of  the  range  resolution  cell  of  interest,  R’  cos  d  =  h.  In  addition,  X^CR-R'jO) 
will  be  approximated  by  its  value  at  the  center  of  the  range  resolution  cell  of  interest, 
Xu(R-R\0)  —  y^(0,0)  =  1.  Furthermore,  it  will  be  assumed  that  RdR  can  be  approximated 
by  R’dR  and  thatcr(R,  0)— cr(R\  0)  in  the  integral.  These  considerations  enable  equation  (5.  15) 
to  be  approximated  by, 

4.k.  -  nR’(AR)  Jct(R’,6)  |  G  (6)  |2kJ"£-  +  ^  sin  w'  cos  (0b  +  0)J  I  d9  (5' 19) 

where  AR  denotes  the  range  resolution  cell  width,  AR  =  c/2  W  . 

Next,  it  is  assumed  that  the  beamwidth  of  G(0)  is  sufficiently  small  so  that  the 
small  angle  approximation  cos  +  0)  ^  cos  0^-0  sin  can  be  used  in  the  integral. 
Because  of  the  periodic  behavior  of  f^(x)  (Fig.  5.  1),  it  follows  that  the  f^  term  in  equation 
(5.  16)  peaks  at  values  of  0  =  0^  given  by, 
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t  p  0,  i  1  f  .  .  • 


(5.20) 


e 

p 


k’  2 vt 
n  A 
2vt 

T“ 


sin  p’  cos  6^  +  p 

sin  (pf  sin  0L 
b 


Furthermore,  because  the  width  of  the  peaks  of  f^(x)  is  ~  ,  the  only  significant 
rangeSof  integration  are  intervals  of  width 

A0  =  ^  s*n  s^n  J  *  (5.21) 

centered  at  the  values  of  0  given  by  equation  (5.  21).  Finally,  it  will  be  assumed  that 
G(0),  f^(x)  and  a  (R\  0)  are  approximately  constant  over  these  integration  intervals. 

Equation  (5.  19)  now  can  be  approximated  by, 


4k.  =  n  R'  (AR)  (A0)  |  il  (0)  |2  p  ?  a  <R\  6p)  |  G  (6  )  |2 


(5.22) 


The  presence  of  the  infinite  sum  makes  equation  (5.  22)  look  quite  formidable; 
however,  in  most  cases  of  interest  this  sum  only  contains  one  significant  term.  As 
indicated  in  section  2,  the  AMTI  system  will  usually  operate  with  its  parameters  chosen 
so  as  to  not  alias  the  clutter.  This  means  that  the  prf  is  chosen  so  that 


h  "  T  $ 


(5.23) 


where  g  £  1  is  a  constant  that  determines  by  how  much  the  clutter  spectrum  is  over  sampled. 
This  means  that, 


2vt  ,  . 

r —  sin  o  sin  9- 
A  b 


ipi 

2vt 


x  sin  6b 


=  I  p  I  *  p 


(5.24) 


where  the  antenna  beamwidth  is  approvimately  given  by  =  X/a  sin  0^. 
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In  other  words,  the  angular  distance  between  any  two  terms  in  the  sum  is  p,  times 
an  integral  number  of  beamwidths.  Since  only  those  terms  within  the  main  lobe  of 
the  beam  contribute  significantly  to  the  sum,  it  follows  that, 

4k -nR'  (AR)(Ae)  I  ^(0)  1 2  cr  (R',  0pI)  I  G(0  ,)  I2 

where  pT  denotes  that  integer  which  makes  0^  closest  to  zero.  The  fact  that  the  clutter 
is  not  aliased  also  results  in  a  further  simplification  because  A0,  as  given  by  equation 
(5.21),  can  be  approximated  by, 

AG  -  —  \£  fi 
n 

which  means  that  A 0  is  on  the  order  of  one  of  a  beamwidth  as  long  as  fi  is  not 
chosen  too  large. 


Making  use  of  these  approximations  in  equation  (5.5)  yields  the  result, 


c  c 

A  similar,  but  somewhat  more  complicated  expression  can  be  obtained  for  the 

case  i1  =  0  but,  for  reasons  to  be  explained  in  the  next  section,  this  case  is  not 

an  interesting  one.  One  final  approximation  will  yield  the  desired  simple  approx- 
k’k1 

imation  to  d.,.,  and  that  is  the  assumption  that  cr  (R',9)  is  approximately  constant 
over  the  entire  beamwidth  of  G(0).  Given  this  assumption,  and  the  approximations 
already  made,  equation  (4.5)  shows  that, 

P  -  E  R’  (AR)  <y  (R\  Q  )  I'  |  G  (9)  |2  dQ 

c  p 

=  ER'  (AR)ct(R',  Qp,)*"1  J*  |  g(6)  |  2  d G 


(5.25) 


(5.26) 


(5.27) 


(5.28) 
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where  the  normalized  gain  function  g(0)  is  defined  by+ 

G(e)  =  \j--1  g  (^_1 0) 

klc * 

The  desired  approximation  to  d  is  now  given  by* 
dkV 

-p—  =  pnc  +Q'  I  fi(°)  I2  s  n  |  g  (^-1  0  )  |2 

where  s  denotes  an  antenna  shape  factor  that  is  independent  of  the  variable  system 
parameters  and  given  by, 

S  =f  Ig  (0)  |2d0 


(5.29) 


(5.30) 


(5.31) 


Actual  calculations  show  that  s  does  not  differ  a  great  deal  from  unity  and  so  plays 
a  relatively  unimportant  part  in  (5.  30). 

2 

The  following  approximate  formula  for  |  v.  |  can  be  obtained  from 

1  K 

equation  (5. 17), 


ViV 


i2a=  |fl(0)  |2 


(5.32) 


where  it  has  been  assumed  that  the  value  of  |  v.f^f  |  for  targets  in  the  exact 
center  of  their  "closest”  spatial -temporal  resolution  cell  is  an  adequate  approximation 
for  the  case  where  a  target  is  off  center.  Substitution  of  equations  (5.  30)  and  (5.  32) 
into  equation  (5.  10)  yields  the  desired  approximate  formula  for  the  improvement 
factor  of  the  DDFT  processor  (valid  for  i'/0) 


I  (v  ,  v  )  ^ 
K  o  ao 


nm  |  fj  (0) 


pnc  +  a  I  V0)  I  sp  I  g(i  eD.)  I 


(5.33) 


+ 


Elementary  antenna  theory  shows  that  g(9)  is  approximately  independent  of 
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If  the  clutter  power  outweighs  the  receiver  noise  power,  equation  (5.33) 


can  be  written  in  factored  form  as  follows, 


I  (v  ,  V  )  Sr  (— )  (- 

o  ao  <y 


sj*  i  g  1  ep.)| 


) 


(5.34) 


These  approximate  formulas  will  be  referred  to  in  the  next  section  in 
order  to  give  insight  into  the  way  that  the  computed  improvement  factors  vary  with 
parameters  such  as  m,n,o'  and  fjt. 
c)  Resolution  and  Ambiguity 

The  resolution  of  a  processor  refers  to  its  ability  to  distinguish  two  desired  targets 
whose  parameters  vq  and  v^  differ  only  by  a  small  amount.  Ambiguity,  on  the  other 
hand,  refers  to  the  processor’s  inability  to  distinguish  two  desired  objects  whose 
parameters  differ  by  a  significant  amount.  (These  ambiguous  parameter  values  are 
sometimes  called  blind  speeds).  For  the  DDFT  processor,  these  characteristics  are 
most  conveniently  studied  by  means  of  the  processor’s  ambiguity  function.  This  function 
is  defined  to  be  the  magnitude  of  the  output  of  a  given  spatial -temporal  resolution  cell 
when  a  desired  target  with  parameters  vq  and  v^  is  present  at  the  input  and  is  denoted 
by  (v()>  vao)where  (i,  k)  denotes  the  indices  of  the  output  resolution  cell  under 
consideration.  This  function  is  seen  to  be  given  by  equation  (5.  16), + 


^ik  <V  ^ao^ 


(5.  35) 


where  it  has  been  assumed  that  t  is  chosen  according  to  equation  (5.  23) 


+The  factor /nm  is  unessential  for  the  present  purpose  and  has  been  dropped. 
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Examination  of  the  sin  mnx/sin  x  term  in  equation  (5.  35)  shows  that  it  has 
peaks  whenever 


v 

o 


X 


v  2d" 


<-s 


+  1),  p=0,  ±  1,  . . . 


Since  the  width  of  the  sin  mx/sin  x  term  is —  ,  it  follows  that  this  term  only 

m 

v  ^ 

contributes  substantially  to  equation  (5,  35)  when  is  within  ±  1/2  of  the 

peaks  given  by  equation  (5.  36).  Similarly  the  f ^(x)  terms  has  peaks  whenever, 


(5.36) 


(—  **  -^)  =  \x  -  (-  +  q),  q=0,  ±  l,  ... 
v  v  a  n 


and  there  is  substantial  contribution  to  equation  (5.  35)  only  when  ^  is 

within  ±  i  —  of  these  peaks. 

2  na  r 

The  resolution  properties  of  the  system  can  be  studied  by  focussing  attention 
on  the  behavior  of  in  the  vicinity  of  one  of  its  peaks.  It  can  be  seen  from  the 
discussion  in  the  last  paragraph  that  X-^  remains  substantially  near  its  peak  value 
as  long  as  the  departures  of  the  parameter  v^  and  v^  -  v^q  from  their  values  at 
the  center  of  the  peak  are  bounded  by  anc^  \  'E5  respectively.  A  region  of 

the  (v^,  v  )  plane  for  which  these  conditions  are  met  is  called  a  spatial -temporal 
resolution  cell  and  has  the  shape  depicted  in  Fig.  5.2. 

Examination  of  Fig.  5.2  reveals  that  two  targets  whose  parameters 
are  equal  to  v  1  v^  and  v^',  v  ,f  respectively  can  appear  in  the  same  spatial - 
temporal  resolution  cell  as  long  as  their  parameters  satisfy 


(5.37) 


and 


I  v  f  -  v  " 
1  o 


I  x 
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v'  -  vM 
ao  ao 
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[it  +^d] 


(5.38) 

(5.39) 
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|  n-i-msi] 


Fig.  5.2.  A  spatial -temporal  resolution  cell. 
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Moreover,  if  either  of  these  conditions  is  violated,  the  two  targets  will  not  appear 
in  the  same  spatial -temporal  resolution  cell  unless  the  differences  between  their 
parameters  becomes  so  large  that  they  lie  in  the  vicinities  of  different  peaks  of 
Xik*  The  quantities  appearing  on  the  right-hand  sides  of  equations  (5.38)  and  (5.39)  are 
called  the  velocity  or  temporal  resolution  and  the  spatial  resolution  respectively.  The 
latter  quantity  can  be  related  to  resolution  in  the  angles  oq  and  Gq  in  an  obvious 
manner. 

In  order  to  study  how  the  resolution  of  the  system  depends  on  the  basic  system 
parameters,  n,  m,  N  and^t ,  it  is  necessary  to  express  the  parameters  d  and  a 
appearing  in  equations  (5.38)  and  (5.39)  in  terms  of  them.  This  is  done  by  means 
of  equations  (2.  22)  and  (2.23)  which  yields  the  results, 


x  _  x  (N-l)  (m-l)+  1  +2n  -  i 


2md  2L 


(N-l)  m  +  1 


(5.40) 


fjA  =  (N-l)  m  +  1  A 
na  n  2md 


A 

2L 


(N-l)  (m-1)  +  1  +  2(j.  - 


m 


(5.41) 


Study  of  the  spatial -temporal  resolution  considered  as  a  function  of  N  and  m  reveals 
that  as  N  increases decreases  steadily  whereas  (■^“j  +  —^ )  decreases  until 
a  critical  value  of  N  is  reached  beyond  which  it  begins  to  increase  again.  The 
critical  value  of  N  is  given  by  N  — and  seems  at  first  glance  a  logical  choice  for  the 
value  of  this  system  parameter.  With  N  set  equal  to^LuT,  it  can  be  seen  that  both  resolu¬ 
tions  degrade  with  increasing  m  thus  leading  to  the  choice  of  m=2.  With  the  velocity  reso¬ 
lutions  optimized  in  this  fashion,  it  turns  out  that  the  resolution  improvement  over  a 
system  with  parameters  N=l,  m=2  is  about  a  factor  of  5.  On  the  other  hand,  the  price 
paid  for  this  resolution  improvement  is  a  reduction  of  Rmax  by  a  factor  of  about  4 
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(for  N=30,  (i=2)  and,  as  will  be  shown  in  tne  next  section,  a  5  db  decrease  in  the  output 
SIR  of  the  system.  In  view  of  the  premium  usually  placed  on  large  values  of  ^max 
and  SIR  in  an  AMTI  system,  it  seems  that  the  increased  resolution  obtainable  by 
making  N  >  1  will  not  be  deemed  worth  the  price  and  that  the  choice  N=l,  m=2  will 
be  the  one  most  often  made. 

A  very  useful  tool  for  visualizing  the  resolution  and  ambiguity  properties  of 
the  processor  is  the  ambiguity  diagram.  This  diagram  indicates  by  means  of 
shading  those  areas  in  the  (v^,  vaQ)  plane  for  which  the  ambiguity  function  y  is 
substantially  different  from  zero.  For  the  case  N=1  and  m  arbitrary  the  discussion 
above  shows  that  the  pertinent  ambiguity  diagram  is  as  indicated  in  Fig.  5.3. 
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Fig.  5.3.  An  ambiguity  diagram  for  N=l. 
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6.  System  Performance 


The  purpose  of  this  section  is  to  present  and  interpret  computed  results  for 
the  MTI  improvement  factor  of  the  multi -antenna  system  under  consideration.  These 
results  were  computed  on  the  basis  of  formulas  derived  in  preceding  sections  descri¬ 
bing  the  signals,  the  interference  statistics,  the  processing  algorithms,  and  the  im¬ 
provement  factor  itself.  As  a  preliminary  to  presenting  the  computed  results,  the 
pertinent  formulas  will  be  summarized  and,  when  necessary,  specialized  to  a  form 
more  suitable  for  computation  or  interpretation.  The  spatial -temporal  factorization 
of  the  improvement  factor,  introduced  in  Section  4d,  will  be  drawn  upon  heavily  for 
interpretation  and  clarification  of  the  results. 

The  first  set  of  results  will  illustrate  the  variation  of  improvement  factor 
with  the  antenna  and  signalling  parameters  introduced  in  Section  2.  These  results 
will  be  drawn  upon  to  suggest  a  reasonable  selection  of  these  parameters.  An  in¬ 
crease  in  the  number  of  antennas  from  one  to  two  will  be  shown  to  produce  dramatic 
improvement,  while  use  of  more  than  two  antennas  provides  little  added 
advantage.  The  effect  of  the  number  of  pulses  transmitted  on  each  antenna,  and  of 
the  parameter  N  governing  the  relationship  between  the  prf  and  the  time  necessary 
to  fly  the  distance  between  phase  centers,  will  also  be  examined.  In  general,  the 
improvement  factor  will  be  plotted  as  a  function  of  radial  target  velocity,  for  a  target 
assumed  located  at  the  antenna  boresight.  This  first  set  of  results  will  be  concluded 
with  a  separate  discussion  of  the  change  of  the  improvement  factor  curves  with  target 
angle  from  boresight. 

The  second  set  of  results  will  illustrate  the  dependence  of  the  improvement 
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factor  on  the  clutter  and  noise  parameters  defined  in  sections  3  and  4.  The  aperture 
effect,  antenna  and  pulse  shape  mismatches,  and  receiver  noise  each  will  be  con¬ 
sidered.  Two  additional  effects,  that  of  slight  inaccuracy  in  the  knowledge  of  air¬ 
craft  velocity,  and  of  internal  motion  of  the  clutter,  will  be  introduced  and  discussed. 

All  the  results  mentioned  above  will  be  presented  for  the  DDFT  processor, 
due  to  the  convenient  implementability  of  this  processor  and  the  fact  that  its  per¬ 
formance  is  very  close  to  the  performance  of  the  optimum  linear  processor.  To 
justify  this  procedure,  computed  comparisons  between  the  DDFT  and  optimum-pro¬ 
cessor  improvement  factors  will  be  presented  at  the  end  of  the  section. 

a)  General  Description  of  Performance  Curves 

The  MTI  improvement  factor  for  the  system  was  represented  in  equation 
(4.  15)  as  9 


where  s  is  the  signal  vector,  w*  is  the  receiver  weighting  vector,  and  M  is  the 
normalized  covariance  matrix  of  the  interference.  This  formula  is  valid  for  any 
linear  processor,  but  will  be  studied  here  only  for  the  DDFT  processor  whose 
weights  are  (ignoring  the  scale  factor  l/^nm  in  (5.  1)) 


w*ik =  WL  exp 


(6.2) 


where  the  weight  w*.^  is  applied  to  the  kth  received  pulse  from  the  ith  antenna,  and 
for  the  optimum  linear  processor  whose  weighting  vector  is 


w  =  M 


s 


(6.3) 
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The  DDFT  processor  is  separable  by  construction.  Therefore  for  both  these  pro¬ 
cessors,  the  argument  of  section  4d  implies  that  subject  to  the  assumption  that  the 
system  performance  is  limited  by  clutter,  rather  than  receiver  noise,  the  improve¬ 
ment  factor  can  be  expressed  as 

I(v  ,  v  )  =  I  (v  )  I  (v  ,  v  )  ,  (6. 4) 

'  o  ao7  s  '  o'  V  o  ao7 

where  the  temporal  improvement  factor  1^.  is  obtained  through  processing  the  pulses 
received  at  each  antenna,  and  the  spatial  improvement  Ig  is  obtained  through  antenna- 
to -antenna  processing. 

The  explicit  dependence  of  I  on  target  radial  velocity  v^  and  on  the  target 
location  parameter  v^Q  is  indicated  in  the  description  of  the  desired  signal  as  (see 
3.  15)). 

s.k  =  exp  (-j2n-2-  i)  exp  [-j2rr(vo  -  v^)  ^  k]  (6.5) 

i  =  0,  1,  . . . ,  m-1 
k  =  0,  1,  . . .  ,  n-  1 

til  til 

where  s.^  represents  the  k  received  pulse  at  the  i  antenna  and  v  depends  on  the 
target  location  angles  (see  (3.  10))  according  to 

v  =  v  sinp  cos  (0  +0,)  .  (6.6) 

ao  vo  v  o  b7  v  7 

Most  of  the  important  features  of  the  system  performance  become  apparent  if  I  is 

plotted  as  a  function  of  v  ,  with  v  =  0.  This  generally  will  be  the  procedure  foll- 

o  a  o 

owed  except  in  section  6.  6.  iv,  where  the  effect  of  variations  of  v  due  to  variations 
K  ao 
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of  target  azimuth  0q  will  be  studied.  The  reason  that  the  spatial  improve¬ 
ment  may  factor  in  (6.4)  depends  only  on  vq  is  apparent  from  (6.5)  where  it  may 
be  noted  that  the  spatial,  or  antenna-to-antenna,  phase  factor  is  independent  of  v  . 
In  the  sequel,  the  notation  I(vq)  will  be  used  for  I(v^,o). 

The  spatial  phase  factor  in  (6.5)  is  periodic  in  v  /v,  with  period  A/2d.  For 

the  case  where  v  =0  and  where  each  antenna  transmits  from  the  same  location, 
ao 

N=1  in  equation  (2.  19)  and  r  =  It  follows  that  the  temporal  phase  factor  has 

period  \/2md  in  v  /v,  and  thus  also  repeats  at  internals  of  X/2d.  These  periodi¬ 
cities  already  have  been  illustrated  in  Fig.  5.3.  In  the  performance  curves  to  be 
presented,  I  will  be  plotted  as  a  function  of  vq/v,  and  usually  only  the  interval 
(0,  \/2d)  will  be  shown.  For  N  =  1,  which  will  be  shown  below  to  be  the  situation 
of  greatest  interest,  I  is  periodic  with  this  period,  and^j-  maybe  straightforwardly 
interpreted  as  the  first  blind  speed  of  the  system.  For  N  ^  1,  the  temporal  phase 
factor  in  (6.5)  has  period  greater  than  X/ 2d ;  this  effect  will  be  discussed  later.  In 
plotting  I  versus  v  /v,  a  grid  consistent  with  the  velocity  sensitivity  of  the  system 
is  obtained  by  changing  vq/v  in  increments  of  ty^mnd),  so  that  mn  points  are  plotted 
in  the  interval  (o,  A/2d).  Reference  to  Fig.  5.3  shows  that  the  velocity  increment 
is  just  sufficient  to  guarantee  the  placing  of  a  target  in  the  center  of  each  DDFT 
resolution  cell.  For  N  ^>1,  analysis  of  Fig.  5.2  shows  that  this  grid  is  still  con¬ 
sistent  with  the  system's  velocity  sensitivity  although  it  could  be  coarsened  without 
losing  any  resolution  in  the  resulting  plots. 

In  an  operating  AMTI  system,  only  those  n(m-l)  resolution  cells  whose 
spatial  indices  i  are  not  equal  to  zero  will  be  examined  for  the  presence  of  moving 
targets.  The  reason  for  this  is  that  a  target  whose  velocity  vq  is  so  small  that  it 
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appears  most  strongly  in  a  resolution  cell  whose  spatial  index  is  zero  cannot  be  dist¬ 
inguished  from  stationary  objects  on  the  ground.  The  following  plots  of  the  improve¬ 
ment  factor  were  obtained  by  first  constructing  the  target  returns  that  result  when 

v  x 

assumes  the  values  -  =  p  ),  p  =  0,  ...  mn-1.  The  improvement  factor 

is  then  computed  for  that  spatial -temporal  resolution  cell  ’’closest”  to  v^  where  the 
term  ’’closest”  is  to  be  taken  in  the  sense  defined  near  the  end  of  Section  5a.  An 
exception  to  this  rule  occurs  when  the  closest  cell  has  indices  (i',k’)  and  i*  =  0;  in 
this  case,  the  improvement  factor  is  computed  for  the  cell  (1,  kf). 

The  improvement  factor  I  depends  not  only  on  the  signal  vector  s  but  also 
on  the  covariance  matrix  M  as  indicated  in  equation  (6.  1).  To  prepare  for  Section 
6c,  where  this  latter  dependence  will  be  examined,  the  structure  of  M  will  be  reviewed 
and  a  specific  form  for  the  antenna  pattern  will  be  introduced  to  allow  for  the  compu¬ 
tation  of  numerical  results.  The  elements  of  M  were  represented  in  equation  (4.  9)  as 


m  ,  =  p  6..,  6  ,,  +  (1  -o'JT.i)  H  (k-k’) 
n  Knc  n  kk*  x  li  v 


(6.7) 


The  function  H(k-k’)  is  the  normalized  correlation  function  of  the  pulses  received 
at  each  antenna.  Now  if  a  Gaussian  illumination  of  each  of  the  antenna  apertures  is 
assumed  and  if  it  is  further  assumed  that  the  range  resolution  is  fine  enough  so  that 
the  elevation  angle  to  the  ground  varies  negligibly  over  a  range  resolution  cell,  it 
can  be  shown  that  the  function  H(k-k’)  has  the  Gaussian  form 

2 


H(k-k')  =  exp  j~2  [£  (k-k')]  j 


(6.8) 
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where 


P=BLL-smo'  (6.9) 

U  ct 

and  where  o'  denotes  the  elevation  angle  from  the  aircraft  to  the  center  of  the  range 
cell.  A  choice  of  the  constant  =  4  yields  a  reasonable  match  between  the  Gaussian 
correlation  function  and  those  which  would  be  obtained  with  other  common  antenna 
illuminations  such  as  uniform  or  triangular  illumination.  The  correlation  function 
(6.8)  is  actually  only  valid  when  the  antenna  beam  is  perpendicular  to  the  aircraft 
velocity,  or  =  90°.  For  6^  ^  90°,  the  right-hand -side  of  (6.8)  must  be  multiplied 
by  a  linear  phase  factor  in  (k-kf)  corresponding  to  the  non-zero  Doppler  frequency 
of  the  clutter  patch  at  the  center  of  the  beam.  However,  this  phase  factor  is  known, 
appears  in  both  the  target  and  clutter  return,  and  can  be  removed  by  suitable  pre¬ 
processing  of  the  data.  For  this  reason,  it  will  be  ignored  in  the  sequel  and  (6.8) 
will  be  used. 

It  may  be  observed  from  (6.  9)  that  /?  is  inversely  proportional  to  the  prf , 
l/r  Thus  a  high  ft  implies  a  prf  which  is  low  in  comparison  to  the  clutter  spectral 
width  and  a  low  fi  corresponds  to  a  higher  prf.  More  detailed  interpretation  of  0 
will  be  given  in  section  6.c.  i. 

b)  Effect  of  Antenna  and  Signalling  Parameters 

Results  now  will  be  presented  which  illustrate  the  effect  of  the  parameters 
m,  n,  and  N  on  the  improvement  factor.  These  results  will  be  presented  for  the 

DDFT,with  a  comparison  between  the  DDFT  and  optimum  processors  being  reserved  for 
section  6.d.  A  target  located  on  the  antenna  boresight  0  =0,  will  generally  be  assumed, 

but  a  separate  discussion  of  the  effect  of  an  off-boresight  target  will  be  given.  It  should 
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be  recalled  from  section  2  that  m  and  N  essentially  determine  d,  a,  t,  and  R  . 

max 

in  studying  the  effects  of  m  and  N,  the  effects  of  changes  in  these  other  parameters 
are  automatically  considered. 

(i)  Number  of  Antennas 

In  Fig.  6.  1,  the  improvement  factor  is  plotted  for  the  case  of  1,  2, 

4  and  8  antennas.  The  parameter  choices  cy  =  10  fi  =  1.0,  =0,  N  =  1  have 

been  made  and  the  product  nm  has  been  fixed  at  t>4.  The  location  parameters  0^ 
and  0Q  have  been  set  equal  to  90°  because  this  choice  does  not  misrepresent  the 
results  as  will  be  demonstrated  in  section  6.c.i.  One  period  of  the  one-antenna 
curve  has  been  plotted  in  the  interval  (0,  A/2d)  for  the  purposes  of  simple  com¬ 
parison  with  the  other  curves,  although  the  parameter  d  is  meaningless  for  the 
one-antenna  case. 

Inspection  of  Fig.  6.  1  shows  that  there  is  a  dramatic  gain  to  be  had 
in  going  from  a  one -antenna  system  to  a  two -antenna  system  in  that  the  velocity 
region  over  which  the  improvement  factor  is  large  enough  to  allow  targets  to  be 
seen  (an  I  of  about  45  db  is  necessary  if  the  input  SIR  is  -30  db)  is  much  larger 
in  the  two  antenna-case.  In  fact  at  the  low  prf  (high  /?)  considered  in  this  example, 
there  is  essentially  no  region  of  visible  velocities  for  the  one-antenna  system. 
Referring  to  equations  (6.  9)  and  (2. 17)  it  is  apparent  that  /?  =  1.0  corresponds  to 
fj.  =2.0.  This  ratio  ^  of  prf  to  nominal  clutter  width  is  not  high  enough  to  produce 
visible  regions  in  the  aliased  clutter  spectrum  seen  at  each  antenna. 

On  the  other  hand,  it  is  apparent  that  even  with  m  =  8,  the  only 
improvement  over  the  m  =  2  case  is  a  slight  narrowing  of  the  blind  regions.  It 
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v,/v 


Fig.  6. 1.  Effect  of  varying  m. 
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seems  that  essentially  all  the  clutter  cancellation  possible  is  achieved  with  2  antennas, 
and  that  little  is  gained  by  increasing  the  number  of  antennas  beyond  2.  The  result 
in  Fig.  6.  1  is  consistent  with  the  approximate  gain  formula,  equation  (5.  34),  which 
indicates  that  I  is  proportional  to  times  a  term  independent  of  the  number  of 
antennas. 

It  has  just  been  demonstrated  that,  as  far  as  the  improvement  factor 
is  concerned,  nothing  substantial  can  be  gained  by  using  more  than  two  antennas  that 
could  not  be  gained  simply  by  increasing  the  number  of  pulses  per  antenna  n.  It  is 
conceivable,  however,  that  the  resolution  of  the  system  could  be  improved  by  in¬ 
creasing  the  number  of  antennas.  Reference  to  Fig.  5.  2  shows  that  the  only  way 
m  influences  the  system  resolution  is  through  the  quantity  md,  which  equation  (2.  22) 
shows  cannot  vary  substantially  from  the  value  L  no  matter  what  choice  of  N  or  m  is 
made.  This  means  that  the  system  resolution  is  not  affected  substantially  and  since 
it  a'ready  has  been  shown  that,  for  N  =  1,  the  improvement  factor  does  not  get  better 
when  m  is  increased,  it  seems  reasonable  to  choose  a  two-antenna  system,  at  least 
when  N  =  1.  It  will  be  pointed  out  in  the  discussion  below  that  a  two -antenna  system 
also  turns  out  to  be  the  most  reasonable  choice  when  N  >  1. 

(ii)  Number  of  Pulses  per  Antenna 

The  effect  of  increasing  n,  while  keeping  the  other  parameters  fixed, 
is  straightforward  and  is  illustrated  in  Fig.  6.  2  for  m  =  2,  N  =  1,  n  =  10  ^ ,  and 
ft  =  1.  Essentially,  doubling  n  adds  3  db  to  the  gain  by  doubling  the  coherent  inte¬ 
gration  time  of  the  processor.  In  principle,  it  seems  desirable  to  make  n  as  large 
as  possible.  In  practice,  n  is  limited  by  such  considerations  as  (a)  when  n  increases 
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Fig.  6.2.  Effect  of  varying  n. 
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acceleration  effects  which  alter  the  phase  structure  of  the  returns  must  be  accounted 
for,  (b)  for  large  n,  the  target  might  have  to  be  tracked  through  several  range  gates 
during  the  coherent  integration  time,  and  (c)  as  the  aircraft  flies  a^ng  with  a  fixed 
beam  pointing  angle,  a  given  target  remains  in  the  beam  only  for  a  finite  period  of 
time. 


(iii)  The  Parameter  N 


The  parameter  N  represents  the  number  of  pulses  transmitted  by  a 


given  antenna  before  the  next  antenna  in  line  moves  to  the  original  spatial  position  of 
that  given  antenna.  It  entered  into  the  clutter  cancellation  constraint  relating  the 
prf  and  the  inter-antenna  spacing  and  aircraft  velocity  according  to  (see  equation 
2.  19) 


VT  =  d  /  (N  -  1+  1/m) 


(6.12) 


Incorporating  (6.  12)  into  (6.  5)  for  the  case  where  v  =0  yields 

ao 


v 


(6.  13) 


The  first,  or  spatial,  phase  factor  in  (6. 13)  is  periodic  in  — -  with  period 


p  .  X 
Ps  _^T 


(6.14) 


while  the  temporal  phase  factor  has  period 


Pt  =  (N  -  1  +  1/m) 


(6. 15) 


V 


It  follows  that  targets  with  — ~  separated  by  an  integer  multiple  of 


Po  =  "23"  [(N-l)  m  +  l] 


(6. 16) 
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V 

are  indistinguishable  and  that  I  repeats  with  this  period  in 

For  N=l,  “"23“*  corresPonding  to  the  repetition  period  of  the  curves  presented 
up  to  now. 

Now  the  effect  of  N  may  be  examined  by  application  of  the  factorization 
result  (6.4),  specialized  to  the  case  v  =  0  as 

10  loSlO  1  (vo>  =  10  loSiO  Js  (vo)  +  10  logio  h  (vo)  (6-  17) 

An  implication  of  (6. 17)  is  that  I  for  an  m-antenna,  n-piPse  system  may  be  obtained 
by  adding  (in  db)  the  improvement  factor  for  an  m-antenna,  1-pulse  system  to  that 
for  a  1-antenna,  n-pulse  system.  To  illustrate  consider  the  example  shown  in  Fig. 

_  c 

6.3.  The  top  trace,  labelled  Ig,  is  the  improvement  factor  for  m  =  2,  n  =  1,  a  10 
Ig  has  period  as  indicated  in  (6*  14).  The  second  curve  corresponds  to  m  =  1, 
n  =  16,  jf?  =  1;  it  has  been  drawn  with  period  ^  ,  which  corresponds  to  N  =  1, 
m  =  2  in  (6.  15)  Adding  the  two  results,  I  is  obtained  for  a  system  where  m  =  2, 
n  =  16,  a  =  10  jS  =  1,  and  N  =  1. 

Figure  6.4  depicts  the  same  situation  as  Fig.  6.3,  except  that  N  =  2.  The 
3  X 

change  is  that  now  so  that  the  horizontal  scale  for  I  has  been  expanded 

with  respect  to  that  for  Ig.  The  overall  period  is  now  Pq  =  3  ,  so  that  to  plot 

one  cvcle  of  the  overall  improvement  factor  curve,  I  must  be  taken  through  3 
cycles  and  1^.  through  2  cycles.  For  N  =  3,  the  results  would  change  in  a  similar 
way.  The  overall  period  would  become  5  ^  ,  corresponding  to  5  cycles  of  I 

and  2  cycles  of  1^.. 
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Fig.  6.  3.  Factorization  of  improvement  factor. 
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Fig.  6.4.  Improvement  factor  for  N=2. 
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Comparing  the  bottom  curves  of  Figs.  b.  3  and  6.  4,  it  is  apparent 

that  increasing  N  produces  essentially  no  improvement  in  system  improvement 

factor.  The  shape  of  the  curve  changes  somewhat,  but  the  peak  and  average 

improvement  factors  stay  about  the  same.  An  explanation  for  this  is  that 

v 

the  variations  of  I  with  — are  essentially  determined  by  the  variations  of  1^, 

since  the  variations  of  I  are  relatively  small.  Thus  the  precise  way  in  which  I 

t  s 

and  1^  line  up  has  an  unimportant  effect  on  the  overall  improvement  factor. 

If  N  were  varied  for  a  fixed  m  >  2,  the  effect  on  the  improvement 

factor  would  be  similar.  The  forms  of  I  and  L  would  stay  the  same  but  the  way 

s  t 

in  which  the  two  line  up  would  change.  However,  there  would  be  no  important 
change  in  the  overall  improvement  factor.  It  m  were  increased  for  a  fixed  N  >  1, 
a  result  similar  to  that  depicted  in  Fig.  6.  1  would  ensue.  The  greatest  improve¬ 
ment  would  be  in  going  from  one  to  two  antennas,  and  the  remaining  improvement 
in  going  from  2  antennas  upward  would  be  small. 

A  loss  not  apparent  in  the  improvement  factor  is  incurred  in  increasing 
N.  As  indicated  in  equation  (2.23),  the  individual  antenna  length,  a,  decreases  with 
N.  This  implies  a  decrease  in  input  signal -to -cl utter  ratio  and  therefore  a  decrease 
in  output  SIR,  even  when  the  improvement  factor  stays  the  same. 

It  was  also  indicated  in  section  2  that  the  unambiguous  range  coverage 
decreases  with  N.  The  only  possible  reason  for  choosing  N  >  1  is  that  slightly  better 
velocity  and  angle  resolution  could  be  achieved,  as  discussed  in  section  5c.  However, 
given  the  fact  that  output  SIR  and  range  coverage  are  of  prime  importance  in  this 
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system,  it  seems  most  reasonable  to  choose  N  =  1  as  will  be  done  in  the  remainder 
of  this  section. 


(iv)  Off -Bore sight  Target 


In  the  results  presented  thus  far,  has  been  taken  to  be  zero,  so 
that  only  the  variation  of  I  with  target  velocity,  and  not  the  variation  with  target 
angle,  has  been  exhibited.  Here  the  effect  of  the  target  angle  from  boresight  0q, 
will  be  considered.  This  effect  is  most  easily  illustrated  for  the  case  of  oq  =  0^  = 
90°,  where  (6.6)  becomes 


v  =  -v  sin  0 
ao  o 


which  may  be  taken  as 


V  wV  0 
ao  o 


(6. 18) 


(6.19) 


since  small  angles  from  boresight  are  not  to  be  considered.  Substituting  (6.  19) 
into  (6.5),  the  signal  return  takes  the  form 

2vt  ,1 


sik=exp 


T  •  Trr  V°  2d  ' 

[-J  21  —  X  1 


exp 


-j  2tt  \  —  +  0 
j  \  V  o 


(6.  20) 


v 

The  temporal  phase  factor,  when  considered  as  a  function  of  — - —  is  simply  advanced 
by  0q.  Thus  a  target  with  zero  velocity,  at  an  angle  0q,  will  give  the  same  temporal 
phase  factor  as  a  boresight  target  with  a  velocity  vq  =  v  0q.  This  angle -velocity 
coupling  effect  is  well-known  in  airborne  doppler  radar  systems,  and  has  been  dis¬ 
cussed  in  more  general  terms  in  section  5c. 
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A  factorization  condition  like  (6.  17)  holds  for  the  off-boresight  case, 
v 

except  that  1^  now  depends  on  both  — - —  and  0q.  In  fact,  since  0^  simply  adds  linearly 
to  the  apparent  velocity  seen  in  the  temporal  phase  factor,  the  counterpart  of  (6.  17) 
is 


10  log10  I  (yn,  ert)  =  10  log,n  IC  (vj  +  10  login  If  (Vn  +V0J 


(6. 17a) 


Here  Ig  is  the  same  as  in  the  boresight  case,  and  1^  (v  +  v0^)  may  be  obtained  by 
simply  shifting  the  result  obtained  in  the  boresight  case. 

Before  presenting  a  result  indicating  the  effect  of  an  off-boresight 

target,  a  warning  mentioned  in  section  4  regarding  the  definition  of  the  improvement 

factor  must  be  recalled.  The  improvement  factor  has  been  defined  as  the  ratio  — 

^in 

of  output  to  input  signal-to-clutter  ratio.  As  is  apparent  from  equation  (4.  1),  the 
input  signal-to-clutter  ratio  decreases  with  increasing  0q,  according  to  the  fall  off 
of  the  antenna  gain  function  with  increasing  0  .  This  causes  output  signal-to-clutter 
ratio  to  generally  decrease  as  the  target  moves  off-boresight,  in  a  manner  which  is 
not  reflected  in  the  behavior  of  the  improvement  factor. 

An  example  of  the  effect  of  0q  on  the  improvement  factor  is  illustrated 
in  Fig.  6.5,  which  corresponds  to  Fig.  6.3  in  all  respects  except  that  an  off-bore¬ 
sight  angle  of  0Q  =  -j  “  has  been  assumed.  Application  of  equations  (2.  18)  and 
(2.  19)  for  the  case  where  —  =  /?  =  1  yields  the  result  a  =  8d,  so  that  I  is  advanced 
by  0q  =  as  indicated.  Comparing  Figs.  6.3  and  6.5,  it  is  observed  that  there 

is  essentially  no  loss  in  improvement  factor,  although  the  shape  of  the  bottom  curves 
are  somewhat  different.  However,  making  the  Gaussian  beam  shape  assumption 
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Fig.  6.5.  Effect  of  off-bore  sight  target. 
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which  led  to  (6.8)  and  (6.  9),  it  is  found  that  for  this  value  of  0  the  input  signal-to- 

o 

clutter  ratio  is  about  5.7  db  down  from  its  value  at  boresight,  and  therefore  the 
output-to-clutter  ratio  is  lower  by  this  same  amount, 
c)  Clutter  and  Noise  Parameters 

The  effect  on  MTI  improvement  factor  of  the  various  clutter  and  noise  para¬ 
meters  which  enter  into  the  interference  covariance  matrix  will  now  be  considered. 
The  interference  effects  to  be  dealt  with  include  (1)  the  aperture  effect,  (2)  antenna 
and  pulse  shape  mismatches,  (3)  receiver  noise,  (4)  inaccurate  knowledge  of  air¬ 
craft  velocity,  and  (5)  internal  motion  of  the  clutter.  The  last  two  of  these  effects 
have  for  simplicity  not  been  included  in  the  interference  matrix  description  developed 
above,  and  will  be  introduced  for  the  first  time  in  this  section.  Again  the  computed 
results  will  be  presented  for  the  DDFT  processor. 

(i)  The  Aperture  Effect 

Reviewing  the  factorization  result  derived  in  section  4  and  restated  in 
(6.4),  it  is  clear  that  the  aperture  effect,  which  influences  only  the  correlation  among 
pulses  seen  at  an  individual  antenna,  affects  only  the  temporal  improvement  factor 
and  not  the  spatial  improvement  factor.  Thus  in  studying  the  aperture  effect,  it 
suffices  to  restrict  attention  to  a  one-antenna  system. 

Recalling  (6.7)  -  (6.9),  and  assuming  again  that  clutter  dominates 
over  white  noise,  the  correlation  coefficient  between  the  kth  and  the  k'th  clutter 
sample  received  by  a  given  antenna  may  be  described  as 
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(b.  21) 


2 


kk 


exp  -j  [fi  (k-k')] 


m  . 


n 


where 


P=P0—  sin<p'^4—  sin  o’ 


(6.22) 


Inspection  of  these  two  equation  can  indicate  qualitatively  the  effect  of  variations  of 
ft  on  the  improvement  factor.  For  very  small  /?,  (6.  21)  indicates  that  the  correlation 
function  of  the  clutter  returns  is  quite  flat.  This  corresponds  to  a  clutter  spectrum 
which  is  narrow  with  respect  to  the  prf,  so  that  very  high  improvement  factor  can  be 
achieved  in  the  frequency  ranges  (clear  regions)  where  the  spectral  amplitude  of  the 
clutter  is  low.  For  higher  /?,  the  clutter  samples  becomes  less  correlated,  the 
clutter  spectral  width  becomes  comparable  to  the  prf,  and  these  clear  regions  in  the 
clutter  spectrum  begin  to  disappear.  The  relationship  between  and  the  prf  is  in¬ 
dicated  explicitly  in  (6.22).  It  may  be  recalled  from  (2.20)  that 


(6.  23) 


corresponding  to  /?  =  2,  specifies  a  prf  equal  to  the  nominal  clutter  spectral  width. 

If  1/t  is  chosen  higher  than  in  (6.23),  /3will  decrease  and  clear  regions  in  the  clutter 
spectrum  will  begin  to  appear.  If  l/r  is  chosen  lower  than  in  (6.23),  the  clutter 
spectrum  will  become  badly  aliased  and  look  essentially  flat. 

In  Fig.  6.6,  the  effect  of  changing  jSfor  a  one-antenna,  16  pulse  system  is 
illustrated.  For  a  one-antenna  system,  the  improvement  factor  has  period  X/2vr 
in  v  /v,  and  one  period  of  this  improvement  factor  is  shown  in  the  figure.  The 
target  is  assumed  to  be  located  at  the  center  of  the  antenna  beam.  Essentially  I(v^) 
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Fig.  6.6.  Effect  of  varying  (3. 
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traces  out  the  inverse  of  the  clutter  spectrum  seen  at  the  antenna.  For  high  ft,  both 
the  clutter  spectrum  and  I(v^)  are  flat.  For  low  ft,  the  clutter  spectrum  becomes 
small  near  a  Doppler  frequency  corresponding  to  half  the  prf  and  the  improvement 
factor  peaks  up  at  this  point.  Two  window  functions  (see  equation  (6.2)),  a  uniform 
window 

w*  =  1  (6. 24) 

and  a  Hamming  window 

=  .54  -  .46  cos  (H*)  (6.25) 

were  tried  for  each  value  of  £  considered  in  Fig.  6.  6.  For  high  ft  where  the  clutter 
spectrum  is  flat,  the  uniform  window  gives  essentially  the  same  performance  as 
the  Hamming  window.  For  low  ft  the  lower  spectral  sidelobes  of  the  Hamming 
window  lead  to  considerably  better  performance  than  with  the  uniform  window  in 
regions  where  the  clutter  spectral  amplitude  is  low.  For  these  cases,  the  results 
plotted  correspond  to  those  obtained  with  the  Hamming  window. 

The  results  plotted  in  Fig.  6.  6  also  illustrate  the  effect  of  changing  beam 
elevation  angle  on  the  MTI  improvement  factor.  As  can  be  seen  from  (6.  22),  the 
effect  of  changing  o’  is  simply  to  change  the  effective  ft  If  the  prf  is  kept  fixed, 
while  o’  decreases  from  90°,  ft  decreases  along  with  the  clutter  spectral  width, 
and  the  effect  of  varying  ft  on  I(v^)  is  as  displayed  in  Fig.  6.  6. 

(ii)  Antenna  Pattern  and  Pulse  Shape  Mismatches 

The  antenna  pattern  mismatches  and  pulse  shape  mismatches  between 
antennas  influence  only  the  correlation  of  clutter  returns  from  antenna-to-antenna, 
and  hence  affect  only  the  spatial  improvement  factor  and  not  the  temporal  improve - 
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ment  factor.  Thus  in  studying  these  effects,  it  suffices  to  consider  an  m-antenna 
system  transmitting  only  one  pulse  per  antenna.  Specializing  (6.7)  to  the  case  of  an 
m-antenna,  one  pulse  per  antenna  system,  where  clutter  dominates  over  white  noise, 
the  correlation  between  clutter  samples  at  the  i  th  and  i’th  antenna  becomes 


m..  =  1  -  S’..,  (6.  26) 

Li  n 

Here  as  defined  in  section  5,  is  a  parameter  measuring  the  antenna -to -antenna 
decorrelation  of  the  clutter  returns  due  to  the  combined  effect  of  antenna  pattern 
and  pulse  shape  mismatches. 

An  analytical  result  for  the  improvement  factor  of  a  two-antenna,  one-pulse 
system  can  be  easily  obtained.  The  covariance  matrix  for  this  case  is 

"l  l-o/“ 


M  = 


Li  -  ck  1  J  . 

This  signal  vector  is  simply  (see  (6.5)). 


(6.27) 


*-C 


/  •  orr  °  2d* 

exp  (-  j  2n  —  -y) 


The  weighting  vector  for  the  DDFT  processor  is 


(6.28) 


(6.29) 


corresponding  to  p  =  1,  m  =  2,  n=  1  in  (6.  2);  the  other  weighting  vector  for  this  system 
(corresponding  to  p  =0)  matches  to  a  zero  velocity  target  and  will  never  be  used  for 
moving  target  detection,  due  to  the  high  correlation  of  the  clutter.  Substituting  these 
last  three  relations  into  the  improvement  factor  formula 
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(6.30) 


KVo)  = 


—% - 

w  M  w 


produces  the  result 

I(v  )  = 
v  o' 


1  ,.>2d  VoN 

1  -  cos  (2n  — - ) 

x  A  v 


<y 


(6.31) 


The  effect  of  rvon  the  improvement  factor  is  readily  apparent  in  (6.31).  The  direct 

relationship  between  the  system  performance  and  the  degree  with  which  antennas 

and  pulse  shapes  can  be  matched  is  clearly  indicated.  In  Fig.  6.7  the  improvement 

-1  _3  -5 

factor  is  plotted  for  m  =  2,  n  =  1,  and  <y  =  10  ,10  ,  and  10  .  Recall  that  for  a 

system  which  transmits  more  than  one  pulse  per  antenna  the  improvement  factor 
would  be  augmented  by  a  simple  addition  (in  db)  of  the  improvement  factor  for  a  one- 
antenna,  n-pulse  system. 

The  effect  of  exon  a  system  with  more  than  two  antennas  is  essentially  the 
same  as  that  described  in  (6.31)  and  Fig.  6.7.  The  proportionality  to  I(vq)  to  l/a 
is  also  valid  for  an  m-antenna,  one-pulse  system.  The  effect  of  increasing  m  for 
a  fixed  o' has  already  been  illustrated  in  Fig.  6.  1. 

(iii)  Effect  of  White  Noise 

The  effect  of  white  noise  is  essentially  to  limit  the  attainable  improve¬ 
ment  factor  of  the  system.  If  clutter  could  be  cancelled  perfectly,  the  MTI  improve¬ 
ment  factor  would  be  limited  by  white  noise.  Once  antenna  patterns  and  pulse  shapes 
have  been  matched  as  closely  as  possible,  the  system  must  be  designed  so  that  the 
receiver  noise  is  low  enough  to  allow  the  improvement  factor  that  the  clutter  cancell¬ 
ation  should  permit.  On  the  other  hand,  given  a  specific  level  of  receiver  noise, 
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Fig.  6.  7.  Effect  of  varying  a. 
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little  is  gained  by  cancelling  clutter  more  precisely  than  to  the  point  where  the  improve¬ 
ment  factor  becomes  limited  by  receiver  noise. 

To  see  this  more  quantitatively,  consider  the  approximate  improvement  factor 
formula  (see  (5.33)) 


nm 


I(v  ,v  )  = 
o  ao 


fi<o)i 


Pnc+«  1  f ,  <0>  1 2  SI*  |  g  (V)  8  ,)|2 


(6.32) 


2 

Here  is  the  input  noise-to-clutter  ratio  per  pulse;  the  factors  | f ^ (0) |  ,  s  and  /i 

- 1  2 

are  all  constants  on  the  order  of  magnitude  of  unity;  and  |g  (>^  )  )  |  ,  the  normalized 

antenna  pattern  function  is  proportional  to  the  clutter  power  present  in  the  output  of  the 
temporal  Doppler  filter  where  the  target  appears. 

The  effect  of  white  noise  will  be  negligible  if 


Pnc  <<(*  |  g  (*_1)  0p.)|2  (6-33) 

and  will  dominate  if 

Pnc»|g(^1)  ep,)|2  .  (6.34) 

In  intermediate  situations  both  the  white  noise  and  the  clutter  affect  the  improvement 
factor.  When  p  is  such  that  the  two  terms  in  the  denominator  of  (6.  32)  are  equal, 
the  white  noise  will  cause  a  3  db  loss  in  improvement  factor. 

In  Fig.  6.8,  the  effect  of  p  is  illustrated.  The  case  considered  is  m  =  2, 
n  =  16,  &  =  10  ,  $  =  1.  For  p  =  4  x  10  ,  the  white  noise  produces  about  a  3  db 

loss  in  improvement  factor  relative  to  the  p  =  0  case  at  the  peaks  of  I(vq).  With 

-7 

p^^  =  10  ,  the  white  noise  has  negligible  effect,  and  the  result  is  indistinguishable 

-4 

from  that  for  p  =0.  For  p  =10  the  output  noise-to-clutter  ratio  is  much 
rnc  rnc  * 
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Fig.  6.8.  Effect  of  white  noise. 
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greater  than  1  and  the  entire  system  is  white  noise  limited. 


(iv)  Imprecise  Knowledge  of  Aircraft  Velocity 
In  order  to  operate  the  multi -antenna  system  in  a  DPCA  mode,  the 
time  offset  between  the  pulse  trains  transmitted  by  adjacent  antennas  must  be  set 
equal  to  d/v  so  that  the  antennas  transmit  from  the  same  points  in  space.  Any  error 
in  knowledge  of  the  aircraft  velocity  will  cause  the  antennas  to  transmit  from  different 
locations  in  space.  This  will  cause  a  decorrelation  between  the  clutter  returns  received 
at  different  antennas,  which  in  turn  will  reduce  the  precision  with  which  clutter  is  can¬ 
celled  and  hence  detract  from  the  system  improvement  factor. 

To  introduce  this  effect  into  the  covariance  matrix,  the  phase  factor 
exp[-  j  2tt|  vaQ  -  A)  i]  (6.35) 

in  (3. 13),  which  was  dropped  in  deriving  (3.  15)  and  in  obtaining  the  description  of 
the  clutter  returns,  must  be  carried  out  through  the  entire  derivation  of  the  clutter 
covariance  matrix.  The  result  is  that  (4.  9)  is  modified  to  become 

mkk,’  =  p  6..,  6  ,  +  (1  -  or  3\.)H  [(k-k*)  +  y  (i-i*)]  (6.36) 

11  rnc  u  kk  '  u  L  1  J 

where  1  , 

Y=:p(A-£)  (6.37) 

is  the  ratio  of  the  timing  error,  A  -  d/v,  to  the  pulse  repetition  time  r.  Assuming 
the  Gaussian  correlation  function  (6.8)  and  neglecting  the  white  noise  yields  the 
equation 

2 

mkk  =  (1  ‘  “S.j)  exP  I"  -5  $2  [(k-k’)  +  Y  (i-i’)]  I  (6.38) 
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The  parameter  y  may  be  related  to  the  percentage  velocity  error  if  it  is  assumed  that 


A  =  d/ (v-Sv) 

where  Sv  is  the  velocity  error.  The  result  is 

Sv  _ 

—  -  my 


(6.39) 


(6.40) 


where  N  =  1  has  been  assumed  so  that  (2.  14)  implies  that  vt  =  md. 

In  Fig.  6.  9  the  effect  of  y  is  illustrated  for  m  =  2,  n=  16,  /9  =  1.0, 

-5  -3 

and  ot  =  10  .  For  y  -  (3.  16)  x  10  a  noticeable  but  small  decrease  of  I  results. 

However,  the  region  of  visible  velocities  is  affected  negligibly.  This  value  of  y  is 

-3 

equivalent  to  a  relative  velocity  error  5v/v  =  (6.32)  10  ,  corresponding  to  a  velocity 

accuracy  readily  achievable  with  present-day  inertial  navigation  systems.  With 
-2 

y  =  10  ,  or  6  v/v  =  .02,  the  visible  velocities  remains  fairly  wide. 

(v)  Internal  Motion  of  the  Clutter 

Internal  motion  of  the  clutter,  such  as  that  caused  by  wind,  will  cause 
a  slight  decorrelation  among  the  clutter  returns  seen  by  the  different  antennas,  and 
thus  limit  the  ability  to  cancel  clutter.  This  is  a  prime  limiting  effect  in  a  ground- 
based  system,  and  an  aim  of  the  DPCA  procedure  is  to  make  the  airborne  system 
perform  as  well  as  a  ground-based  system.  Only  when  the  airborne  system  is 
achieving  this  objective,  and  thus  performing  quite  well,  will  internal  clutter  motion 
be  of  importance. 

To  avoid  additional  complexity,  the  analysis  leading  to  the  clutter 
covariance  matrix  (4.  9)  assumed  stationary  scatterers  on  the  ground.  If  the  clutter 
scatterers  are  assumed  to  have  a  distribution  density  A(R,  0,v  )  in  velocity  as  well 


90 


Fig.  6.  9.  Effect  of  unknown  aircraft  velocity. 


91 


as  in  range  and  azimuth  and  if  an  independence  assumption  analogous  to  equation  (3.  21) 
is  made,  then  the  resulting  normalized  covariance  matrix  (neglecting  receiver  noise) 
is  given  by 

m^,  =  (1-0-6. .,)  H  (k-k’)  Lj^p  (k-k')  +  (i-i') }  (6.41) 

where  OJ  v 

'c> 


L(x)  =■ 


JCT(vc)exp|-j2n^  xj 


dv 


(6.42) 


J  <r(v  )  dv 


c  c 

In  arriving  at  equation  (6.41)  it  has  been  assumed  that  the  clutter  cross  section 
density  cr(R,  0,  v  )  can  be  factored  in  the  form  cr(R»  0,  v  )  =  cr(R,  0)  cr  (v^).  It  also 
should  be  noted  that  the  total  clutter  input  power  per  pulse  now  is  given  by 
equation  (4.  5)  multiplied  by  J*  cr(vc)  dv^.  Equation  (6.41)  shows  that  internal  clutter 
motion  enters  as  a  simple  multiplicative  effect  in  the  covariance  matrix. 

In  order  to  deal  with  (6.41)  numerically,  the  correlation  function 
L(x)  will  be  assumed  to  be  composed  of  a  dc  component  due  to  the  zero-velocity 
scatterers  plus  a  Gaussian  component  due  to  moving  scatterers,  so  that 

L(x)  -  Kj  +  K2  exp  (-  \  (ex)2  )  .  (6. 43) 

The  effect  of  internal  clutter  motion  is  illustrated  in  Fig.  6. 10.  For 
-  2  -2 

K2/K1  =  10  and  e  =  4.  46  x  10  about  a  3  db  decrease  in  improvement  factor  results. 
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Fig.  6. 10.  Internal  motion  of  clutter. 
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d)  Comparison  of  DDFT  and  Optimum  Linear  Processor 


A  series  of  examples  will  now  be  presented  to  compare  the  improvement 
factors  achieved  with  the  DDFT  and  the  optimum  linear  processors.  It  was  shown 
theoretically  in  Section  5  that  the  DDFT  processor  should  perform  nearly  as  well  as 
the  optimum,  and  the  computed  results  bear  this  out.  In  practice,  the  DDFT  pro¬ 
cessor  would  be  used  because  of  its  simpler  implementation. 

In  Fig.  6.  11,  the  two  processors  are  compared  for  m=l,  n=lb,  and  0=1, 
and  the  results  are  essentially  indistinguishable.  As  in  the  one-antenna  curve  of 
Fig.  6.  1  I(vq)  is  relatively  flat  because  the  prf  has  been  chosen  low  enough  so  that 
the  aliased  clutter  spectrum  is  relatively  flat.  With  a  flat  clutter  spectrum,  little 
is  gained  through  the  precise  tailoring  of  the  Doppler  filter  sidelobes  achievable 
with  the  optimum  processor.  Fig.  6.  12  compares  the  two  processors  for  m=l,  n=lb  and 
4.  A  Hamming  window  was  used  with  the  DDFT  processor.  Here  the  prf  has  been 
chosen  at  five  times  the  nominal  clutter  width  (see  (6.  10)  and  (6. 11))  so  that  at  Doppler 
frequencies  near  one -half  the  prf  the  clutter  spectral  amplitude  is  quite  low  and  the 
improvement  factor  is  high.  In  this  frequency  region  the  optimum  processor  does 
produce  substantially  higher  improvement  factor  than  the  DDFT  processor.  But  this 
difference  is  considered  unimportant  for  two  reasons.  First,  in  the  frequency  range 
where  the  difference  is  appreciable  the  DDFT  processor  is  already  performing  so 
well  that  targets  will  be  detected  and  additional  improvement  factor  is  useless. 

Secondly,  it  is  desirable  in  this  AMTI  system  to  use  as  low  a  prf  as  possible  to 
maximize  range  coverage;  hence,  the  situation  depicted  in  Fig.  6.  11  is  more  likely 
than  that  in  Fig.  6.  12. 
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IMPROVEMENT  FACTOR  (dB)  ^  IMPROVEMENT  FACTOR  (dB) 


ig.  6.  11.  Comparison  of  DDFT  and  optimum  processors. 


Fig.  6. 12.  Comparison  of  DDFT  and  optimum  processors. 
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In  Fig.  6.  13  the  two  processors  are  compared  for  m=2,  n=l,  ot  =  10 
The  results  are  identical,  and,  in  fact,  it  can  be  shown  theoretically  that  for  this 
case  the  DDFT  processor  is  identical  to  the  optimum  linear  processor.  This  impor¬ 
tant  result  may  be  generalized  to  state  that  in  an  m-antenna  system  where  ot  and 
are  the  only  parameters  of  importance  in  the  interference  covariance  matrix,  the 
spatial  improvement  factors  for  the  DDFT  and  optimum  linear  processors  are 
identical. 

In  Fig.  6.  14,  m=2,  n=16,  ^=10  /3=1  and  p  =  4  x  10  so  that  the  effect 

nc 

of  white  noise  is  included.  Obviously  the  DDFT  processor  produces  results  essentially 
equivalent  to  the  optimum  processor. 

In  Fig.  6.15,  the  effect  of  all  interference  parameters  is  included.  Again  the 
performance  of  the  two  processors  is  indistinguishable. 

7.  Conclusion 

The  preceding  sections  have  developed  a  basic  theory  of  multi-antenna  AMTI 
radars  and  used  this  theory  to  compare  various  system  configurations.  The  main  con¬ 
clusion  of  this  work  is  that,  subject  to  certain  assumptions  a  two-antenna  system  with 
each  antenna  emitting  alternate  equally  spaced  pulses  (N=l)  offers  a  great  improvement 
over  a  conventional  single -antenna  radar. 

The  main  assumption  on  which  this  conclusion  is  based  is  that  antenna  mis¬ 
matches  and  pulse  shape  mismatches  can  be  controlled  well  enough  so  that  values  of 

-4  -5 

the  degradation  parameter  ot  on  the  order  of  10  -  10  can  be  achieved.  Future  work 

should  be  directed,  in  part,  to  establishing  a  better  mathematical  and  physical  under¬ 
standing  of  the  factors  that  contribute  to  ot.  Another  topic  for  further  investigation 
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ig.  6.  13. 


Comparison  of  DDFT  and  optimum  processors. 


Fig.  6.  14.  Comparison  of  DDFT  and  optimum  processors. 
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is  the  comparison  of  system  configurations  on  the  basis  of  constant  illuminated  range, 


R  ,  rather  than  on  the 
max 

Preliminary  work  on  such 
arrived  at  on  a  constant  p 
R  basis. 


max 


basis  of  a  constant  sampling  rate,  jxt  as  was  done  above, 
a  comparison  seems  to  indicate  that  the  general  conclusions 
basis  also  hold  when  comparisons  are  made  on  a  constant 


Fig.  6.  15.  Comparison  of  DDFT  and  optimum  processors. 
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